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Whereas semiclassical gravity is based on the semiclassical Einstein equation with sources given by 
the expectation value of the stress-energy tensor of quantum fields, stochastic semiclassical gravity 
is based on the Einstein-Langevin equation, which has in addition sources due to the noise kernel. 
The noise kernel is the vacuum expectation value of the (operator-valued) stress-energy bi-tensor 
which describes the fluctuations of quantum matter fields in curved spacetimes. We show how 
a consistent stochastic semiclassical theory can be formulated as a perturbative generalization of 
semiclassical gravity which describes the back reaction of the lowest order stress-energy fluctuations. 
The original approach used in the early investigations leading to the establishment of this fleld 
Q (for a review emphasizing ideas, see [^) was based on quantum open system concepts (where the 
metric field acts as the "system" of interest and the matter fields as part of its "environment") and 
the infiuence functional method. Here, following Refs. we first give an axiomatic derivation 
of the Einstein-Langevin equations and then show how they can also be derived by the original 
method based on the influence functional. As a first application we solve these equations following 
Ref. and compute the two-point correlation functions for the linearized Einstein tensor and 
for the metric perturbations in a Minkowski background. We then turn to the important issue of 
the validity of semiclassical gravity by examining the criteria based on the ratio of the variance of 
fluctuations of the stress-energy tensor of a quantum field to its mean. We show a calculation of 
these quantities performed in Ref. for a massless scalar field in the Minkowski and the Casimir 
vacua as a function of an intrinsic scale defined by introducing a smeared field or by point separation. 
Contrary to prior claims, the ratio of variance to mean-squared being of the order unity does not 
necessarily imply the failure of semiclassical gravity. Expressions for the variance to mean-squared 
ratio as a function of the ratio of the intrinsic to extrinsic scale (defined by the separation of the 
plates or the periodicity of space for the Casimir topology) identifying the spatial extent where 
negative energy density prevails are useful for studying quantum field effects in worm holes, baby 
universe and the design feasibility of 'time-machines'. This study also raises new questions into the 
very meaning of regularization of the stress-energy tensor for quantum fields in curved or ordinary 
flat spacetimes. It corroborates the view held by one of us that at high energies it is more useful to 
consider spacetime as possesing an extended structure, with the stress-energy bi-tensor of quantum 
matter fields as the source. To compare with the more familiar point-defined quantum field theory 
operative at low energies, one needs to address the issue of regularization of the noise kernel. This 
we discuss following the work of Ref. which uses the point-separation method to derive a general 
expression for a regularized noise-kernel for quantum fields in general curved spacetimes. From 
these expressions describing the behavior of fiuctuations of quantum fields one could investigate a 
host of important problems in early universe and black hole physics. This includes addressing the 
viability of a vacuum-dominated phase which infiationary cosmology is predicated upon, examining 
how structures are seeded, how the horizon behavior of a black hole is altered and how the late time 
dynamics of black holes is affected with backreaction of Hawking radiation. We intend to discuss 
these applications in a later report. 
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I. OVERVIEW 



Stochastic semiclassical gravity (STG) is a theory developed in the Nineties using semiclassical gravity (SCG, 
quantum fields in classical spacetimes, solved self-consistently) as the starting point and aiming at a theory of quantum 
gravity (QG) as the goal. While semiclassical gravity is based on the semiclassical Einstein equation with the source 
given by the expectation value of the stress-energy tensor of quantum fields, stochastic gravity (we will often use 
the shortened term stochastic gravity as there is no confusion in the source of stochasticity in gravity being due to 
quantum fields here, and not of a classical origin) includes also its fluctuations in a new stochastic semiclassical or 
the Einstein-Langevin equation. If the centerpiece in semiclassical gravity theory is the vacuum expectation value 
of the stress-energy tensor of a quantum field, and the central issues being how well the vacuum is defined and how 
the divergences can be controlled by regularization and renormalization, the centerpiece in stochastic semiclassical 
gravity theory is the stress-energy bi-tensor and its expectation value known as the noise kernel. The mathematical 
properties of this quantity and its physical content in relation to the behavior of fluctuations of quantum fields 
in curved spacetimes are the central issues of this new theory. How they induce metric fluctuations and seed the 
structures of the universe, how they affect the black hole horizons and the backreaction of Hawking radiance in black 
hole dynamics, including implications on trans-Planckian physics, are new horizons to explore. On the theoretical 
issues, stochastic gravity is the necessary foundation to investigate the validity of semiclassical gravity and the viability 
of inflationary cosmology based on the appearance and sustenance of a vacuum energy-dominated phase. It is also 
a useful beachhead supported by well-established low energy (sub-Planckian) physics to explore the connection with 
high energy (Planckian) physics in the realm of quantum gravity. 

In these lectures we want to summarize our work on the theory aspects since 1998. (A review of ideas leading to 
stochastic gravity and further developments originating from it can be found in Ref. [y^]; a comprehensive formal 
description is given in Refs. It is in the nature of a progress report rather than a review. In fact we will have 

room only to discuss two or three theoretical topics of basic importance. We will try to mention all related work so 
the reader can at least trace out the parallel and sequential developments. The references at the end of each topic 
below are representative work where one can seek out further treatments. 

Stochastic gravity theory is built on three pillars: general relativity (GR), quantum fields (QF) and nonequilibrium 
(NEq) statistical mechanics. The first two uphold semiclassical gravity, the last two span statistical field theory, or 
NEqQF. Strictly speaking one can understand a great deal without appealing to statistical mechanics, and we will try 
to do so here. But concepts such as quantum open systems || l O | and techniques such as the influence functional [Tl| ] 
(which is related to the closed-time-path effective action |l2| , |l3| ) were a great help in our understanding of the physical 
meaning of issues involved towards the construction of this new theory, foremost because quantum fluctuations and 
correlation have become the focus. Quantum statistical field theory and the statistical mechanics of quantum field 
theory [p"4|-pT[ also aided us in searching for the connection with quantum gravity through the retrieval of correlations 
and coherence. We show the scope of stochastic gravity as follows: 



I. Ingredients 

A. From General Relativity ||l^,|l9| to Semiclassical Gravity. 

B. Quantum Field Theory in Curved Spacetimes pC|-[2^: 
Stress-energy tensor: Regularization and renormalization. 
Self-consistent solution: Backreaction problems ||2^ . 

Effective action: Closed time path, initial value formulation [T^JT^ . 
Equation of motion: Real and causal. 

C. Nonequilibrium Statistical Mechanics: 
Open quantum systems 

Influence Functional: Stochastic equations | pl| . 

Noise and Decoherence: Quantum to classical transition |p5| , p6| . 

D. Decoherence in Quantum Cosmology and Emergence of Classical Spacetimes [^Tj 

II. Theory 

A. Dissipation from Particle Creation p8| , ^ . 
Backreaction as Fluctuation-Dissipation Relation (FDR) | |30t . 

B. Noise from Fluctuations of Quantum Fields [P j3^j3l|] . 
* C. Einstein-Langevin Equations 



2 



* D. Metric Fluctuations in Minkowski spacetime 



III. Issues 



* A. Validity of Semiclassical Gravity 
B. Viability of Vacuum Dominance and Inflationary Cosmology. 

* C. Stress- Energy Bitensor and Noise Kernel: Regularization Reassessed 

IV. Applications: Early Universe and Black Holes 

[Test field treatment:] 

A. Wave Propagation in Stochastic Geometry |3^ . 

B. Black Hole Horizon Fluctuations: Spontaneous/ Active versus Induced/Passive ^4 
[Including backreaction:] 

C. Noise induced inflatio n ||39tl - 

D. Structure Formation: [|I0|, Trace anomaly-driven inflation [^-Q. 

E. Black Hole Backreaction as FDR 



V. Related Topics 

A. Metic Fluctuations and Trans- Planckian Problem [|^,|9|j50). 

B. Spacetime Foam [ p]| . 

C. Universal 'Metric Conductance' Fluctuations [ |5^ . 

VI. Ideas 



A. General Relativity as Geometro- Hydrodynamics [ p3[ . 

B. Semiclassical Gravity as Mesoscopic Physics [Q. 

C. From Stochastic to Quantum Gravity: 

via Correlation hierarchy of interacting quantum flelds j|,|l^,|l^ . 
Possible relation to string theory and matrix theory . 

We list only the latest work in the respective topics above describing ongoing research. The reader should consult 
the references therein for earlier work and the background material. We do not seek a complete coverage here, but 
will discuss only the selected topics in theory and issues (Parts II and III above), i.e., those topics marked with *s 
above, in Sections 2-5 below. In a later report we plan to discuss recent advances in the applications (Parts IV-V 
above). We use the (+, +, +) sign conventions of Refs. [|l^,^, and units in which c = h = 1. 



II. STOCHASTIC SEMICLASSICAL GRAVITY 



In this section we introduce stochastic semiclassical gravity. We will start with semiclasical gravity which assumes 
that the stress-energy fluctuations of the matter flelds are negligible. When these fluctuations cannot be neglected 
we may perturbatively correct the semiclassical theory to take them into account. This will lead to the semiclassical 
Einstein-Langevin equations ^ as the dynamical equations which describe the back reaction of quantum stress-energy 
fluctuations on the gravitational fleld. We will introduce these equations by two independent methods. The flrst, 
based on Ref . [Q , is axiomatic and motivated by the search of consistent equations to describe the back reaction of the 
quantum stress-energy fluctuations on the gravitational fleld. The second, based on Ref. [^, is by functional methods. 
It is motivated by the observation [|l| that in some open quantum systems classicalization and decoherence [p5[ | on the 
system may be brought about by interaction with an environment. The environment being in this case the matter 
fields and some degrees of freedom of the quantum gravitational field [p^ . 



A. Semiclassical gravity 



Semiclassical gravity describes the interaccion of the gravitational field assumed to be a classical field with quantum 
matter fields. This theory cannot be derived as a limit in a quantum theory of gravity, but can be formaly derived in 
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several ways. One of them is the leading 1/iV approximation of quatum gravity interacting with N independent and 
identical free quantum fields |]56|] . When one assumes that the fields interact with gravity only and keep the value 
of NG finite, where G is Newton's gravitational constant, one arrives at a theory in which formaly the gravitational 
field can be treated as a c-number (i.e. quatized at tree level) but matter fields are fully quatized. 

In the semiclassical theory theory the expectation value of the stress-energy tensor in a given quantum state is the 
source of the gravitational field. This theory may be summarized as follows. Let (A^,g^i/) be a globaly hyperbolic 
four-dimensional spacetimc with metric 5^1/ and consider a real scalar quantum field cf) of mass m propagating on it 
(we assume a scalar field for simplicity but this is not essential). 

The classical action Sm for this matter field is given by the functional 

where is the covariant derivative associated of the metric 5^^, ^ is a couplig parameter coupling the field to the 
scalar curvature R, and g under the square root stands for the determinat of g^,j. 

The field may be quantized in the manifold using the standard canonical formalism [po|-^ . The field operator in 
the Heisenberg reperesentation is an operator valued distribution solution of the Klein-Gordon equation, the field 
equation of Eq. (|2.lD, 

(□ - - e^)0 = 0. (2.2) 

We will write the field operator as (l)[g]{x) to indicate that it is a functional of the metric. 

The classical stress-energy tensor is obtained by functional derivation of this action in the usual way T'^'^{x) = 
{2/y^)SSjn/Sgfj,u, leading to 

V,,V, + G^,)02, (2.3) 

where □ — V^V^ and G^jy is the Einstein tensor. With the notation T^i,[g,(/)] we explicity indicate that the stress- 
energy tensor is a functional of the metric gp^u and the field 0. 

The next step is to define a stress-energy tensor operator TpAg]{x). Naively one would replace the classical field 
<j) in the above functional by the quantum operator (l)[g\, but this procedure involves taking the product of two 
distributions at the same spacetime point. This is ill-defined and we need a regularization procedure. There are 
several regularization methods which one may use, one is the point-splitting or point-separation regularizaton method 
in which one introduces a point ?/ in a neighborhood of the point x and then uses as a regulator the vector tangent 
at the point x of the geodesic joining x and y; this method is discussed in detail in section ^ Another well known 
method is dimesional regularization in which one works in arbitray n dimensions, where n is not necessarily an integer. 



and then uses as the regulator the parameter e = n — 4; this method is used in this section and in section III. The 



regularized stress-energy operator using the Weyl ordering prescription, i.e. symmetrical ordering, can be written as 

= , V>[g]} + P^''[5]0^[5], (2.4) 

where V'^'^lg] is the differential operator 

2?r ^ - i) 9^''(^)°- + e (R'^i^) - V^V^) . (2.5) 

Note that if dimensional regularization is used, the field operator propagates in a n-dimensional spacetime. Once 
the regularization prescription has been introduced a renormalized and regularized stress-energy operator may be 
defined as 

T^. [9] (x) = [9] (x) + F^, [g] {x)i, (2.6) 

where / is the identity operator and -^^[5] are some symmetric tensor counterterms which depend on the regulator 
and are local functionals of the metric. Here we are assuming that all these terms depend on the regulator, see Ref. ||] 
for details. These states can be chosen in such a way that for any pair of physically acceptable states, i.e. Hadamard 
states in the sense of Ref. [Q, and \(p) the matrix element {i^lT^^ip) defined as the fimit of the previous expression 
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when the regulator takes the physical value, is finite and satisfies Wald's axioms |2|,^. These counterterms can be 
extracted from the singular part of a Schwinger-DeWitt series [^,^ 59 1. The choice of these counterterms is not 
unique but this ambiguity can be absorbed into the renormalized coupling constants which appear in the equations 
of motion for the gravitational field. 

The semiclassical Einstein equations for the metric g^i, can then be written as 



inG 



{G^4g]+Ag^,) 



(2.7) 



where {T^^)[g] is the expectation value of the operator (2.6) after the regulator takes the physical value in some 
physically acceptable state of the field on {A4,gf^^). Note that both the stress tensor and the quantum state are 
functionals of the metric, hence the notation. The parameters G, A, a and f3 are the renormalized coupling constants, 
respectively, the gravitational constant, the cosmological constant and two dimensionless coupling constants which 
are zero in the classical Einstein equations. These constants must be understoood as the result of "dressing" the 
bare constants which appear in the classical action before renormalization. The values of these constants must be 
determined by experiment. The left hand side of Eq. (pj) may be derived from the gravitational action 
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167rG 



(R - 2 A) + aCf^^p^G^^'P^ + PR 



(2.8) 



where C^^pcr is the Weyl tensor. The tensors A^^ and B^^ ^ come from the functional derivatives with respect to the 
metric of the terms quadratic in the curvature in Eq. (2.8), they are explicitly given by 



A,. 



B„ 



1 



(2.9) 
(2.10) 



are the Riemann and Ricci tensors, respectively. These two tensors are, like the Einstein and 
metric tensors, symmetric and divergenceless: V^A^j/ = = V^B^^. Note that a classical stress-energy can also be 



where R^j,vpa and R^^ 



added to the right hand side of Eq. (2.7), but for simplicity we omit such a term. 

A solution of semiclassical gravity consists of a spacetime (Al, 5^1^), a quantum field operator (l)[g\ which satisfies 
( p.l[ ), and a physically acceptable state \'4')[g] for this field, such that Eq. (2.7) is satisfied when the expectation value 
of the renormalized stress-energy operator is evaluated in this state. 

For a free quantum field this theory is robust in the sense that it is consistent and fairly well understood. As long 
as the gravitational field is assumed to be described by a classical metric, the above semiclassical Einstein equations 
seems to be the only plausible dynamical equation for this metric: the metric couples to mater fields via the stress- 
energy tensor and for a given quantum state the only physically observable c-number stress energy-tensor that one 
can construct is the above renormalized expectation value. However, lacking a full quantum gravity theory the scope 
and limits of the theory are not so well understood. It is assumed that the semiclassical theory should break down at 
Planck scales, which is when simple order of magnitude estimates suggest that the quantum effects of gravity should 
not be ignored because the energy of a quantum fluctuation in a Planck size region, as determined by the Heisenberg 
uncertainty principle, is comparable to the gravitational energy of that fluctuation. 

The theory should also break down when the fluctuations of the stress-energy tensor are large. This has been 
emphasized by Ford and colaborators [ |60| , |6ll . It is less clear, however, how to quantify what a large fluctuation is, 
and different criteria have been proposed |6l|,|62PJ^. In Sec. IV we will discuss at length the issue of the validity of 



the semiclassical theory. One may illustrate the problem by the following example inspired in Ref. ||60| . 

Let us assume a quantum state formed by an isolated sytem which consists of a superposition with equal amplitude 
of one conflguration with mass Mi and another with mass Af2- The semiclassical theory as described by Eq. (2.7) 
predicts that the gravitational field of the system is produced by the averaged mass (Afi -f M2)/2. However, one 
would expect that if we send a succession of test particles to probe the gravitational field of the above system half of 
the time they would react to the field of a mass Mi and half of the time to the field of a mass M2. If the two masses 
differ substantially the two predictions are clearly different, note that the fluctuation in mass of the quantum state 
is of order of (Mi — M2)^. Altough this example is suggestive a word of caution should be said in order not to take 
it too literaly. In fact, if the previous masses are macroscopic the quantum system decoheres very quickly [ |63[ and 
instead of being described by a pure quantum state it is described by a density matrix which diagonalizes in a certain 
pointer basis. For observables associated to such a pointer basis the matrix density description is equivalnet to that 
provided by a statistical ensemble. In any case, though, the results will differ from the semiclassical prediction. 
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B. Axiomatic route to stochastic semiclassical gravity 



1. The noise kernel 

The purpose of stochastical semiclassical gravity is to go beyond the semiclassical theory and account for the 
fluctuations of the stress-energy operator. But first, we have to give a physical observable that describes these 
fluctuations. To lowest order, these fluctuations are obviously described by the following bi-tensor, constructed with 
the two-point correlation of the stress-energy operator, 



4N^.pA^,y) = 77({V(a:),tpa(2/)})b], (2.11) 



where the curly brackets mean anticommutator, and where 

if,,{x) = f^,{x)- {f^,{x)). (2.12) 

This bi-tensor will be called noise kernel from now on. Note that we have defined it in terms of the unrenormalizcd 
stress-tensor operator Tfj_^[g]{x) on a given background metric g^^, thus a regulator is implicitly assumed on the r.h.s. 



of Eq. (2.11). However, for a linear quantum field the above kernel is free of ultraviolet divergencies because the 



ultraviolet behaviour of (T^^ (a;)Tpcr(y)) is the same as that of {Tfj_v{x)){Tp„{y)), thus following Eq. (2.6) one can 



replace by the renormalized operator in Eq. ( 2.11 ) ; an alternative proof of this result is given in section 0, 



The noise kernel should be thought of as a distribution function, the limit of coincidence points has meaning only in 
the sense of distributions. An analysis of the noise kernel based on the point-separation method is given in section 
The bi-tensor N^^p„{x,y) is real and positive-semidefinite, as a consequence of T"^ being self-adjoint. A simple 

proof can be given as follows. Let be a given quantum state and let Q be a selfadjoint operator, = 
then one can write { i^\QQ \tp) ~ (V'lQ^QI^) = > 0- Now let i{x) be a selfadjoint operator representing the 



operator t^^ of Eq. (2.11) but where the coordinate x now carries also the tensorial indices of Eq. (2.11), then if 
we define Q — J dxf{x)i{x) for an arbitrary well behaved "function" /(x) (which carries also tensorial indices), the 
previous unequality can be written as J dxdyf{x){^p\t{x)t{y)\ip)f{y) > 0, which is the condition for the noise kernel 
to be posite semi-definite. Note that when we will consider, later on, the inverse kernel N~J-p^{x,y) we will implicitly 
assume that we arc working in the subspace obtained from the eigenvectors which have strictly positive eigenvalues 
when we diagonalize the noise kernel. 

2. The Einstein- Lang evin equation 
Our purpose now is to perturbatively modify the semiclassical t heor y. Thus we will assume that the background 



spacetime metric f;^^ is a solution of the semiclassical Einstein Eqs. (2.7). The stress-energy tensor will generally have 
quantum fiuctuations on that background spacetime. One would expect that these fluctuations will have some effect 
on the gravitational field and that this may now be described by g^^ + hf^^^, where we will assume that ft-^i^ is a linear 
perturbation to the background solution. The renormalized stress-energy operator and the state of the quatum field 
may now be denoted by [g + h] and \tlj)[g + h], respectively, and (T^) [g + h] will be the corresponding expectation 
value. 

Let us now introduce a Gaussian stochastic tensor field defined by the following correlators: 

(^^-(a;)). = 0, {^pAx)^pAy))s = Nf,,p^{x,y), (2.13) 

where (. . .)s means statistical average. The symmetry and positive semi-definite property of the noise kernel guarantees 
that the stochastic field tensor i^^^ just introduced is well defined. Note that this stochastic tensor does not capture 
the whole quantum nature of the fiuctuations of the stress-energy operator since it assumes that cumulants of higher 
order are zero, but it does so at the Gaussian level. 

An important property of this stochastic tensor is that it is covariantly conserved in the background spacetime 
V'^^^ijy = 0. In fact, as a co nsequence of the conservation of T/^j^[g] one can seen that '^'^N p^^pa{x,y) = 0. Taking 



the divergence in Eq. ( |2.13D one can then show that (V^^^,/)s = and (V^^^,y(a::)V(J^pCT(y))s = so that V^^^u is 
detreministic and represents with certainty the zero vector field in M. 

One can also see that for a conformal field, i.e. a field whose classical action is conformally invariant, is traceless, 
i.e. g^'^^fj.p = 0, so that for a conformal matter field the stochastic source gives no correction to the trace anomaly. 
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In fact, from the trace anomaly result which states that g^^Tj^^[g\ is in this case a local c-immber functional of g^^ 



times the identity operator, we have that g'^'^ {x)Nfj^i,p„{x,y) = 0. It then follows from Eq. ( p.l3[ ) that {g^" £,^j.v) s = 



and {g^" {x)£,f_iv{x)g'"^ {y)£,pa(y)) s = 0; an alternative proof based on the point-separation method is given in section 0. 

All these properties make it quite natural to incorporate into the Einstein equations the stress-energy fluctuations 
by using the stochastic tensor as the source of the metric perturbations. Thus we will write the following equation. 



SttG 



(G^^b + /i] + A(.g^^ + h^,^)) - 2iaAf,^ + (3B^^)[g + h] 



{T^.)[g 



(2.14) 



This equation known as the semiclassical Einstein- Langevin equation, is a dynamical equation for the metric perturba- 
tion h^i, to linear order. It describes the back-reaction of the metric to the quantum fluctuations of the stress-energy 
tensor of matter fields, and gives a first order correction to semiclassical gravity as described by the semiclassical Ein- 
stein equation ( ^.7[ ). Note that the stochastic source ^^j/ is not dynamical, it is independent of /i^^ since it decribes 
the fluctuations of the stress tensor on the semiclassical background gfj.^. 

An important property of the Einstcin-Langevin equation is that it is gauge invariant. In the sense that if we change 
hf_ii, h ' + V u.rji, + V uTjp where 77^ is a stochastic vector field on the manifold M. All the tensor which appear in 
the Eq. (2.14) transform as R^i,[g + h'] — Rfj,iy[g + h] + £riRfj.i^[g] t o line ar order in the pert urbations, where is the 
Lie derivative with respect to 77'^. If we substitute h by h' in Eq. ( |2.14| ), we get Eq. ( 2.14 ) plus the Lie derivative of 
the combination of the tensor which appear in Eq. (2.7). This last combination vanishes when Eq. (2.7) is satisfied, 



as we have assumed. It is thus necessary that the background metric g^i, be a solution of semiclassical gravity. 



A solution of Eq. ( p. 14 ) can be formally written as h^^[S]. This solution is characterized by the whole family of 



its correlation functions. From the statistical average of this equation we have that g^^, + {hp,^) g must be a solution 
of the semiclassical Einstein equation linearized around the background g^v The fluctuation of the metric around 
this average are described by the moments of the stochastic field h^f^v[^] — h^u[£,] — {h^,y)s. Thus the solutions of the 
Einstein-Langevin equation will provide the two point metric correlators {hf^j^{x)h^p^{y)) s ■ 

The stochastic theory may be understood as an intermediate step between the semiclassical theory and the full 
quantum theory. In the sense that whereas the semiclassical theory depends on the point-like value of the sress-energy 
operator, the stochastic theory carries information also on the two-point correlation of the stress-energy operator. 

We should also emphasize that, even if the metric fluctuations are classical (stochastic), their origin is quantum not 
only because they are induced by the fluctuations of quantum matter, but also because they are supposed to describe 
some remnants of the quantum gravity fluctuations after some mechanism for decoherence and classicalization of the 
metric field p^-|68|. From the formal assumption that such a mechanism is the Gell-Mann and Hartle mechanism 
of environment-induced decoherence of suitably coarse-grained system variables |64 6^, one may, in fact, derive the 
stochastic semiclassical theory Nevertheless, that derivation is of course formal, given that, due to the lack of 
the full quantum theory of gravity, the classicalization mechanism for the gravitational field is not understood. 



3. A toy model 



To illustrate the relation between the quantum and stochastic descriptions it is useful to introduce the following 
toy model. Let us assume that the gravitational equations are described by a linear field h coupled to a scalar source 
T[(/)^] independent of h, with the classical equations D/i = T in fiat spacetime. We should emphasize that this model 
would describe a linear theory of gravity analogous to electromagnetism and does not mimic the linearized theory of 
gravity in which T is also linear in h. However it captures some of the key features of linearized gravity. 

In the Heisenberg representation the quantum field h satisfies 

Uh = f. (2.15) 

The solutions of this equation, i.e. the field operator at the point x, hx, may be written in terms of the retarded 
propagator Gxy as, 

K^hl + J dx'Gxx'fx', (2.16) 

where is the free field. From this solution we may compute, for instance, the following two point quantum correlation 
function (the anticommutator) 

\{{K,hy}) = \{{hi:hl}) + J dx'dy'Gxx'Gyy,{{fx-,fy,}), (2.17) 
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where the expectation value is taken in the quantum state in which both fields (j) and h have been quantized and we 
have used that for the free field = 0. 

Let us now compare with the stochastic theory for this problem. As discussed previously we should write a classical 
equation for the field h where on the r.h.s. we substitute T by the expectation value (T) plus a Gaussian stochastic field 
^ defined by the quantum correlation function of T (the noise kernel) . Note that in our model since T is independent 
of h we may simply renormalize the expectation value u sing t ime ordering, then for the vacuum state of the field (p, 
we would simply have (f)o = 0. Thus ^ is defined as in (|2a|) by (^), = and {^^Qs = (1/2) ({f^, f^}) - {f^)(fy). 
The semiclassical stochastic equation is thus 



ah = (T)+C 

The solution of this equation may be written in terms of the retarded propagator as, 



dx'G^ 



(2.18) 



(2.19) 



from where the two point correlation function for the classical field h, after using the definition of ^ and that {h^)s = 0, 
is given by 



(2.20) 



Comparing (2.17) with (2.20) we see that the respective second terms on the right hand side are identical provided 
the expectation values are computed in the same quantum state fo r the field (f>, note that we have assumed that T 
does not depend on h. The fact that the field h is also quantized in ( 2.17 ) does not change the previous statement. In 
the real theory of gravity T, in fact, depends also on h and then the previous statement is only true appro ximately, i.e 
perturbatively in h. The nature of the first terms on the right hand sides of equations ( ^.17 ) and ( 2.20 ) is different: 
in the first case it is the two point expectation value of the free quatum field Hq whereas in the second case it is the 
average of the two point classical average of the homogeneous field /iq, which depends on the initial conditions. Now 
we can still make these terms to be equal to each other if we assume for the homogeneous field h a distribution of 
initial conditions such that {h'^hy)s — (l /2)( |/t°, hy}). Thus, under this assumption on initial conditions for the field 
h the two point correlation function of ( 2.20| ) equal the quantum expectation value of ( 2.17 ) exactly. 



Thus in a linear theory as in the model just described one may just use the statistical description given by (2.1 



to compute the quantum two point function of equation (2.16). This does not mean that we can recover all quantum 
correlation functions with the stochastic description, see Ref. for a general discussion about this point. Note 
that, for instance, the commutator of the classical stochastic field h is obviously zero, but the commutator of the 
quantum field h is not zero for timelike separated points; this is the prize we pay for the introduction of the classical 
field ^ to describe the quantum fluctuations. Furthermore, the statistical description is not able to account for the 
graviton-graviton effects which go beyond the linear approximation in h. 



C. Functional approach 



The Einstein-Langevin equation (2.14) may also be derived by a method based on functional techniques Func- 
tional techniques have a long history in semiclassical gravity. It started when effective action methods, which are 
so familiar in quantum fleld theory, were used to study the back-reaction of quantum fields in cosmological models 
Jnl . These methods were of great help in the study of cosmological anisotropics since they allowed the introduction 
of familiar perturbative treatments into the subject. The most common method, the in-out effective action method, 
however, led to equations of motion which were not real because they were taylored to compute transition elements 
of quantum operators rather than expectation values. Fortunately the appropriate technique had already been de- 
veloped by Schwinger and Keldysh |l^ in the so called Closed Time Path (CTP) or in-in effective action. These 
techniques were soon adapted to the gravitational context and were applied to different problems in cosmology. 
As a result one was now able to deduce the semiclassical Einstein equations by the CTP functional method: starting 
with an action for the interaction of gravity with matter fields, treating the matter fields as quantum fields and the 
gravitational field at tree level only. 

Furthermore, in this case the CTP functional formalism turns out to be related | |3l| , |7^ -[75|,pt to the influence 
functional formalism of Feynman and Vernon ]7^ ] since the full quantum system may be understood as consisting 
of a distinguished subsystem (the "system" of interest) interacting with an environment (the remaining degrees of 
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freedom). The integration of the environment variables in a CTP path int egral y ields the influence functional, from 
which one can define an effective action for the dynamics of the system 30 T^JtS 



In our case, we consider the metric field g^v{x) as the "system" degrees of freedom, and the scalar field and 
also some "high- momentum" gravitational modes [|68| as the "environment" variables. Unfortunately, since the form 
of a complete quantum theory of gravity interacting with matter is unknown, we do not know what these "high- 
momentum" gravitational modes are. Such a fundamental quantum theory might not even be a field theory, in which 
case the metric and scalar fields would not be fundamental objects Thus, in this case, we cannot attempt to 
evaluate the influence action of Feynman and Vernon starting from the fundamental quantum theory and performing 
the path integrations in the environment variables. Instead, we introduce the influence action for an effective quantum 
field theory of gravity and matter [8C|j30|] , in which such "high-momentum" gravitational modes are assumed to have 
been already "integrated out." Adopting the usual procedure of effective field theories one has to take the 

effective action for the metric and the scalar field of the most general local form compatible with general covariance: 

S[g, (f)] = Sg[g\ + Sr,^[g, 0] + • • • , 



where Sg[g\ and Sm[g, (f\ are given by Eqs. (2.8) and respectively, and the dots stand for terms of order higher 

than two in the curvature and in the number of derivatives of the scalar field. Here, we shall neglect the higher 
order terms as well as self-interaction terms for the scalar field. The second order terms are necessary to renormalize 
one-loop ultraviolet divergencies of the scalar field stress-energy tensor, as we have already seen. Since is a 
globally hyperbolic manifold, we can foliate it by a family of i = constant Cauchy hypersurfaces S*. We denote by 
X the coordinates on each of these hypersurfaces, and by ti and tf some initial and final times, respectively. The 
integration domain for the action terms must be understood as a compact region hi of the manifold , bounded by 
the hypersurfaces St^ and Et^. 



Assuming the form (2.1) for the effective action which couples the scalar and the metric fields, we can now introduce 
the corresponding influence functional. This is a functional of two copies of the metric field that we denote by g~^i, 
and g~jj. Let us assume that, in the quantum effective theory, the state of the full system (the scalar and the metric 
fields) in the Schrodinger picture at the initial time t = ti can be described by a factorizable density operator, i.e., 
a density operator which can be written as the tensor product of two operators on the Hilbert spaces of the metric 
and of the scalar field. Let p^iti) be the density operator describing the initial state of the scalar field. If we consider 



the theory of a scalar field quantized in the classical background spacetime {M.,gfiu) through the action (2.1), a state 
in the Heisenberg representation described by a density operator p[g] corresponds to this state. Let {[(^(xjp} be the 
basis of eigenstates of the scalar field operator 0'^(x) in the Schrodinger representation: 0^(x) |(/?)^ = </'(^) WY ■ The 
matrix elements of p^{ti) in this basis will be written as pi\i~p, (p\ = ^{^\ P^(ti) 1'^)^- We can now introduce the influence 
functional as the following path integral over two copies of the scalar field: 

TiAg^^g-] ^ 1 p.[</>+(t.),0-(iO] 5[<P+{tf)~ct>-{tf)] (2.21) 



The above double path integral can be rewritten as a closed time path (CTP) integral, namely, as an integral over 
a single copy of field paths with two different time branches, one going forward in time from ti to t/, and the other 
going backward in time from tf to ti. From this influence functional, the influence action of Feynman and Vernon, 
SYp[g^ ,g~], is defined by 

TiY[g+,g-] = e^^'^^o*^a~\ (2.22) 

this action has all the relevant information on the matter fields. Then we can define the effective action for the 
gavitational field, Scfi[g^ ,g^\, as 

S,s[g+,g-] EE Sg[g+] - Sg[g-] + 5if[.9+, .g"]- (2.23) 

This is the action for the classical gravitational fleld in the CTP formalism. However, since the gravitational flcld is 
only treated a t tree level, this is also the effective "classical" action from where the classical equations can be derived. 



Expression (2.21) is ill-deflned, it must be regularized to get a meaningful influence functional. We shall assume that 



we can use dimensional regularization, that is, that we can give sense to Eq. (2.21) by dim ensional continuation of all 
the quantities that appear in this expression. For this we substitute the action Sm in (2.21) by some generalization to 
n spacetime dimensions, and similarly for [5]. In this last case the parameters G, A, a and /3 of Eq. ( |2.8| ) are the 
bare parameters G_b, A^, ub and /3b, and instead of the square of the Weyl tensor one must use {2/?>){R^up„R^^''P" — 
R^iyR^^) which by the Gauss-Bonnet theorem leads the same equations of motion as the action ( p.8[ ) when n = 4. 
The form of Sg in n dimensions is suggested by the Schwinger-DeWitt analysis of the ultraviolet divergencies in the 
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matter stress-energy tensor using dime nsional regularization. One can then write the effective action of Feynman 
and Vernon, Sos[g~^,g~], in Eq. ( 2.23 ) in dimensional regularization. Since both Sm and Sg contain second order 
derivatives of the metric, one should also add some boundary terms ||l^,^. The effect of these terms is to cancel out 
the boundary terms which appear when taking variations of Scff[g~^ , g^] keeping the value of g'^^^ and fixed on the 
boundary oiU. Alternatively, in order to obtain the equations of motion for the metric in the semiclassical regime, 
we can work with the action terms without boundary terms and neglect all boundary terms when taking variations 
with respect to g^^. From now on, all the functional derivatives with respect to the metric will be understood in this 



sense. 



Now we can derive the semiclassical Einstein equations (2.7). Using the de finiti on of the s tress-energy tensor 



g)SSm/Sgfii, and the definition of the influence functional, Eqs. (2.21) and (2.22), we see that 



SSiF[g^,g~ 



\/-g{x) ht'^ix) 



(2.24) 



where the expectation value is taken in the n-dimensional spacetime generalization of the state described by p[g\. 
Therefore, differentiating Scslg'^ ■, g~] in Eq. (2.22) with respect to 17+^, and then setting = g~j^ = g^^, we get the 
semiclassical Einstein equation in n dimensions. This equation is then renormalized by absorbing the divergencies 
in the regularized ( T^'^ ) [g] into the bare parameters, and taking the limit n ^ 00 we get the physical semiclassical 
Einstein equations (|2.7|). 



D. Influence functional route to stochastic semiclassical gravity 



In this subsection we derive the semiclassical Einstein-Langevin equation (2.14) by means of the influence functional. 
We also work out the semiclassical Einstein-Langevin equations more explicitly, in a form more suitable for specific 
calculations. 

In the spirit of the previous derivation of the Einstein-Langevin equa tion s, we now seek a dynamical equation for 
a linear perturbation h^n to a semiclassical metric g^i^, solution of Eq. (2/7). Strictly speaking if we use dimensional 
regularization we must consider the n-dimensional version of that equation; see Ref. |^ for details. From the result of 
the previous subsection, if such equation were simply a linearized semiclassical Einstein equation, it could be obtained 
from an expansion of the effective action S(.s[g + h'^ ,g + h^]. In particular, since, from Eq. ( p^ ), we have that 



{fr{x))[g + h] = 



SSiF [g + h'^, g+h 



y^det(g+7i)(a?) ShJ^uix) 



(2.25) 



h+=h-=h 



the expansion of {TJ^") [g + h] to linear order in h^^i, can be obtained from an expansion of the influence action 
Swig + ,g + ^ ] up to second order in /i^^^. 

To perform the expansion of the influence action, we have to compute the flrst and second ord er fun ctional derivatives 
of SiF[g^ , g~] and then set g'^i, = gj^i, = g fj.,^ ■ If we do so using the path integral representation ( 2.21 ), we can interpret 
these derivatives as expectation values of operators. The relevant second order derivatives are 



1 



s'sMg+,g- 



\/-g{x)^/-g{y) 5gtu{x)5gt^{y) 
1 5^Sw[g+,g-] 



y/-g{x)^y~g{y) 5gZv{x)5gpa{.y) 



-Hr'^[g]{x,y) - K^^''P-[g]{x,y)-V^N^^''P-[g]{x,y), 
~H^''[g]{x.y)-^N^'"'"[g]{x,y), (2.26) 



where 



N^'^P-'[g]{x,y) ^ i {{i^%x), iP%y)]) [g], Hr'^[g]{x,y) ^ \ Im (T*(f '^-(x)f "-(y)) 



Hr"\g\{x,y) ^ -1 ([fM^(x), fP^{y)\)\g\ K^^P-^\g\{x,y) 



^~g{x)^-g[y) \Sg^Ax)Sgpaiy) 



S'S„,[g,^] 



)[g], 



with t'^" defined in Eq. (2.12), and where we use a Weyl ordering prescription for the operators in the last of these 
expressions. Here, [ , ] means the commutator, and we use the symbol T* to denote that, first, we have to time order 
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the field operators (p and then apply the derivative operators which appear in each term of the product T^'^{x)TP'^{y), 
where T'^" is the functional (^j). This T* "time ordering" arises because we have path integrals containing products 
of derivatives of the field, which can be expressed as derivatives of the path integrals whi ch do not contain such 
derivatives. Notice, from their definitions, that all the kernels which appear in expressions ( 2.26 ) are real and also 
that H^'^'"^ is f ree o f ultraviolet divergencies in the limit n 4. 

From Eqs. ( ^.26 ) it is clear that the imaginary part of the influence action, which does not contribute to the 
semiclassical Einstein equation (2.7) because the expectation value of T^^'^[g] is real, contains information on the 
fluctuations of this operator. From (2.24) and ( 2.26| ), taking into account that SiF[g,g] = and that S'if[5~, g"*"] = 
—Sjp[g^,g^], we can write the expansion for the influence action SiF[g + h'^ , g+h'^] around a background metric gat 
in terms of the previous kernels. Taking into account that these kernels satisfy the symmetry relations 

Hr''^ix,y) = Hr^%y,x), Hr^i^.v) = ~HP''^%y,x), K^^^^x^y) = KP^^%y,x), 

and introducing the new kernel 

H^^^P%x,y) ^ Hr'"'ix,y)+Hr'"'ix,y), 
the expansion of Sif can be finally written as 

1 



(2.27) 
(2.28) 



SM9+h+,g + h-] = - I d^x ^^g(x) {T^'%x))[g] [h^,{x)] 



1 

~2 
I 

+ 2 

where we have used the notation 



d'xd^y^^)./^) [V(^)] {H^'''"'[g]{x,y) + K^^^P"[g]{x,y)){hp„{y)} 



d^xd^y^^^)^^) [h^Ax)]N^^''P''[g]{x,y) [V(2/)] + 0(/i'), 



(2.29) 



(2.30) 



We are now in the position to carry out the formal derivation of the semiclassical Einstein-Langevin equation. The 
procedure is well known pl| , ^0|j7^j8^ , it consists of deriving a new "stochastic" effective action using the the following 
identity: 

^_ 1 J d^xd^v^-g(x)^-g{y) [h^^(x)] N ^^-^ {x ,y) [h.^iy)] ^ J pj^j .pj^j ^2 J d^x (x) [h^^x)] ^ ^2.31) 

where V [^] is the probability distribution functional of a Gaussian stochastic tensor ^'^'^ characterized by the correlators 
(2.13) with iV"'"^'' given by Eq. (2.11), and where the path integration measure is assumed to be a scalar under 
diffeomorphisms of (A^,g^y). The above identity follows from the identification of the right hand side of ( p. 31 ) with 
the characteristic functional for the stochastic field In fact, by differentiation of this expression with respect to 
[hjjjX it can be checked that this is the characteristic functional of a stochastic field characterized by the correlators 
( 2.13 ). When N'^'^p"' {x,y) is strictly positive definite, the probability distribution functional for ^^"^ is explicitly given 
by 



^-iJd^xdSy/^^)y/^^e-{^)N-^\A^,y)i'''{y) 



(2.32) 



where N^J-p^[g]{x,y) is the inverse of N^^''P'^[g]{x,y) defined by 

d'zV^)N^-^-ix,z)N:^^\Az,y) = \ (<5,^C + 5^<5,^) ^!£=£l. 

We may now introduce the stochastic effective action as 

Sl^ [h+,h-;g;i] = Sg[g + h+]-Sg[g + h-]^ 5fp [h+ ,h- ; g;(,], 
where the "stochastic" influence action is defined as 



(2.33) 



(2.34) 
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Sfp[h+,h-;g;Cn] = Rc SM9 + h+ , g+h'] + x ^/^^ ^^^^ (x) [h^,{x)] + 0{h^) 



(2.35) 

Note that the influence functional as defined from the influence action ( 2.29| ) can be written as a statistical average 
over ^f^'^: 



(2.36) 



Thus, the effect of the imaginary part of the influence action ( 2.29 ) on the corresponding influence functional is 
equivalent to the averaged effect of the stochastic source ^^'^ coupled linearly to the perturbations h^^^ 
The stochastic equations of motion for h^i, can now be derived as 

1 SS'^f,[h+,h-;g;^, 



y/- dct{g + h)(x) Sh+,ix) 



(2.37) 



h+=h-=h 



Then, from ( [2.25 ), taking into account that only the real part of the influence action contributes to th e exp ectation 
value of the stress-energy tensor, we get, to linear order in /i^j^ the stochastic semiclassical equations ( ^.14 ). To be 
precise we get first the regularized n-dimensional equations with the bare parameters, and where instead of the tensor 
Af"" we get (2/3)^'' where 



(Px- 



'9 {RparujR^' 



1 



_ 



(2.38) 



Of course, when n = A these tensors are related, A^^'^' = (2/3)1)'^''. After that we renormalize and take the limit n ^ A 
to obtain the Einstein-Langevin equations in the physical spacetime. 



1. Explicit linear form of the Eintein-Langevin equation 

We can write the Einstein-Langevin equation in a m ore e xplicit form by working out the expansion of {T'^'^)[g+h] 
up to linear order in the perturbation h^^. From Eq. (2.25), we see that this expansion can be easily obtained from 
( ^.291 ). The result is 

{fr{x))[g + h] = (fr(a;))[g] + {f^'^''''[g;h]{x))[g] ~ 2 [d-y V^)H^'''"^[g]{x,y)hpAy) + Oih'), (2.39) 



where the operator j'^i)'^" is defined from the term of first order in the expansion of T°''[g + h] as 

T'^'ig+h] = T'"'[g] + T^'^^'^ig; h] + 0{h^), (2.40) 

using, as always, a Weyl ordering prescription for the operators in the last definition. Here we use a subscript n on 
a given tensor to indicate that we are explicitly working in n-dimensions, as we use dimensional regularization, and 
we also use the superindex (i) to generally indicate that the tensor is the first order correction, linear in /i^^, in a 
perturbative expansion around the background g^^. 

Using the Klein-Gordon equation ( ^.2| ), and expressions (2.3) for the stress-energy tensor and the corresponding 
operator operator, we can write 

n'^'^ig; h] = Q g'^-hp, - 5>;k - <5;/i(^) rr[<?] + t^^W, h] il[g], (2.41) 

where J-'^'^[g; h] is the differential operator 

= (^£.-2^(h'"'~l9'"''hP^^ + I [VV/i^ + VV^/i^ - ah'"' - VVh^ - g'^VVkp^ 

+ g'^'ahP + [Vh^p + Vhfp ~ Wph"" - 2g'"V''hp^ + g'"V^l)VP - g'^'hp^V'V] . (2.42) 
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It is understood that indices are raised with the background inverse metric g'^'^ and that all the covariant derivatives 

are associated t o the metric g^^,. 

Substituting ( ^.3£ ) into the n-dimensional ver sion of the Einstein-Langevin Eq. ( 2.15 ), taking into account that 
g^jy satisfies the semiclassical Einstein equation (2.7), and substituting expression ( 2.41 ) we can write the Einstein- 
Langevin equation in dimensional regularization as 



SttGb 



g'^'CP" hp„ + G^'PK + + Kb ( h^"' - - g"" hP 



(x) 



- as l^D^'^P" - - gf^DP^hp^ + Df^Ph^p + D^Phi^j {x) - 2pB l^B''^"''-- gP^BP^hp^ + B^Ph^p + B^Ph^ ) {x 
^^2 ^ 2 / ry [g] {x, y) hp, (y) = 2/.- (-4)^'^ i^) , 



(2.43) 

where the tensors G^'', D^''^ and B^'^ are computed from the semiclassical metric g^^, and where we have omitted 
the functional dependence on g^^ and hp^ in G*^"^", D'-'-^p'^, j^wt"^ a^d TP'^ to simplify the notation. The parameter 
/i is a mass scale which relates the dimensions of the phys ical field with the dimensions of the corresponding field 
in n-dimensions, 4>n = Notice that, in Eq. ( 2.43| ), all the ultraviolet divergencies in the limit n ^ 4, which 

must be removed by renormalization of the coupling constants, are in (0^(x)) and the symmetric part RP^p'^ {x,y) 
of the kernel H^'^p'^ {x,y), whereas the kernels N^'^p'^ {x,y) and RP^^p'^ {x^y) are free of ultraviolet divergencies. If we 
introduce the bi-tensor F^'^ p'^ [g]{x , y) defined by 



FP'^P^[g]{x,y)^{iP%x)iP:{y))[g] 
where i^'^ is defined by Eq. ( ^.12 ), then the kernels N and Ha can be written as 



(2.44) 



NP-P" [g] (a;, y) = - Re FP^p- [g] (x, y) , 



y) - 4 Im F^P'^ [g] {x, y), 



(2.45) 



where we have used that 2 {tP" (x) tP'^ (y)) = {{tP^ix), iP'^{y)}) + {[tP^ix), fP'^iy)]), and the fact that the first 
term on the right hand side of th is ide ntity is real, whereas the second one is pure imaginary. Once we perform the 
renormalization procedure in Eq. ( 2.43| ), setting n — A will yield the physical semiclassical Einstein-Langevin equation. 
Due to the presence of the kernel RP'^p"' {x,y), this equation will be usually non-local in the metric perturbation. 



2. The kernels for the vacuum state 

When the expectation values in the Einstein-Langevin equation are taken in a vacuum state |0), such as, for instance, 
an "in" vacuum, we can go further. Since then we can write these expectation values in terms of the Wightman and 
Feynman functions, defined as 

G+{x,y) EE {0\Mx)My)m9], «a„(x,2/) = (0| T(0„(a:)0„(y)) |0)[<?]. (2.46) 

These expressions for the kernels in the Einstein-Langevin equation will be very useful for explicit calculations. To 
simplify the notation, we omit the functional dependence on the semiclassical metric gp^,, which will be understood 
in all the expressions below. 



From Eqs. (2.45), we see that the kernels NP'^p'^{x, y) and RP^'^p'^ {x, y) are the real and imaginary parts, respectively, 
of the bi-tensor FP'^p" {x, y). From the expression ( ^.4[ ) we see that the stress-energy operator TP'^ can be written as 
a sum of terms of the form |yla;^„(a;), Bx4>n{x)^ , where Ax and Bx are some differential operators. It then follows 
that we can express the bi-tensor FP'^p"'{x, y) in terms of the Wightman function as 

Fr'"'{x,y) = VPVPGt{x,y)yy;Gtix,y)+Vy;Gt{x,y)yyPyGt{x,y) 

-|-2PrKG+(x,y)V-G+(x,y)) +2P^-(V^G+(a;,y)V^G+(:r,y)) +2Pr2?r(Gn'(^,y)), (2-47) 



where VP'^ is the differential operator (^.5|). From this expression and the relations ( 2.45| ), we get expressions for the 
kernels Nn and i?A„ in terms of the Wightman function G^{x,y). 

Similarly the kernel RPJ^P'^{x, y), can be written in terms of the Feynman function as 
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- g^%x)V^PG,Sx, y)\^lG,^{x, y) - g''-(2/)VJ!VJG,„(x, 2;)V^X-^a„(.T, y) 
+ \ 9^-ix)g'"'iy)V:V^yG,^ix, y)V,X''a„(x, y) + /Cr(2V^a„(a;, y)V^a„(a;, y) 

- g'^'^ixMG.Jx, y)V:G,S^, y)) + 2 K.^JK.p'^ {G^x, y)) 
where JC'^'^ is the differential operator 

/Cr ^ C (5'^'^(x)n. - V^^ + G^'^ix)) ~ \ m^g^^{x). 



(2.48) 



(2.49) 



Note that, in the vacuum state |0), the term ^)\{x)) in equation ( 2.43| ) can also be written as (0^(a;)) = iGp^{x,x) = 

Finally, the causality of the Einstein-Langevin equation ( 2.43| ) can be e xplici tly seen as follows. The non-local terms 
in that equation are due to the kernel H{x,y) which is defined in Eq. (2.28) as the sum of Hs{x,y) and HA{x,y). 
Now, when the points x and y are spacelike separated, (jinix) and 4>n{y) commute and, thus, G^{x,y) = iGpJ^x,y) = 
(1/2)(0| {(j)n{x) , (j)n{y)} |0), which is real. Hence, from the above expressions, we have that Hj^^'"'{x, y) = Hj;l^'"^{x, y) = 
0, and thus H}j;'^P'^{x, y) = This fact is not surprising since, from the causality of the expectation value of the stress- 
energy o perat or, we know that the non-local dependence on the metric perturbation in the Einstein-Langevin equation, 
see Eq. (2.14), must be causal. 



3. Discussion 

In this section, based on Refs. |^,|| we have shown how a consistent stochastic semiclassical theory of gravity can be 
formulated. This theory is a perturbative generalization of semiclassical gravity which describes the back reaction of 
the lowest order stress-energy fluctuations of quantum matter fields on the gravitational field through the semiclassical 
Einstein-Langevin equation. We have shown that this equation can be formally derived with a method based on the 
influence functional of Feynman and Vernon, where one considers the metric field as the "system" of interest and the 
matter fields as part of its "environment" . An explicit linear form of the Einstein-Langevin equations has been given 
in terms of some kernels which depend on the Wightman and Feynman functions when a vacuum state is considered. 

In Ref. the fluctuation-dissipation relations and particle creation in the stochastic gravity context was also 
discussed. When the background solution of semiclassical gravity consists of a stationary spacetime and a scalar field 
in a thermal equilibrium state, it is possible to identify a dissipation kernel in the Einstein-Langevin equation which 
is related to the noise kernel by a fluctuation-dissipation relation ||83| , p^ . The same relation was previously derived by 
Mottola p6| using a linear response theory approach, and it is in agreement with previous findings on a Minkowski 



background |[7^; see, however, comments in |47|. This result was also generalized to the case of a conformal scalar 
field in a conformally stationary background solution of semiclassical gravity. 

Particle creation by stochastic metric perturbations in stationary and conformally stationary background solutions 
of semiclassical gravity was also considered. It is possible to express the total probability of particle creation and the 
number of created particles (the expectation value of the number operator for "out" particles in the "in" vacuum) in 
terms of the vacuum noise kernel. Remarkably the averaged value of those quantities is enhanced, over the semiclassical 
result, by the presence of stochastic metric fluctuations. In the particular cases of a Minkowski background and a 
conformal field in a spatially fiat RW background, the energy of the created particles can be expressed in terms of 



the vacuum dissipation kernels; see also Refs. [pi 85 



III. METRIC FLUCTUATIONS IN MINKOWSKI SPACETIME 



In this section we describe the first application of the full stochastic semiclassical theory of gravity, where we 
evaluate the stochastic gravitational fiuctuations in a Minkowski background. In order to do so, we first use the 
method developed in section || to derive the semiclassical Einstein-Langevin equation around a class of solutions 
of semiclassical gravity consisting of Minkowski spacetime and a linear real scalar field in its vacuum state, which 
may be considered the ground state of semiclassical gravity. Although the Minkowski vacuum is an eigenstate of the 
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total four-momentum operator of a field in Minkowski spacetime, it is not an eigenstate of the stress-energy operator. 
Hence, even for these solutions of semiclassical gravity, for which the expectation value of the stress-energy operator 
can always be chosen to be zero, the fluctuations of this operator are non- vanishing. This fact leads to consider the 
stochastic corrections to these solutions described by the semiclassical Einstein-Langevin equation. 

We then solve the Einstein-Langevin equation for the linearized Einstein tensor and compute the associated two- 
point correlation functions. Even though, in this case, we expect to have negligibly small values for these correlation 
functions for points separated by lengths larger than the Planck length, there are several reasons why it is worth 
carrying out this calculation. 

On the one hand, these are the first back-reaction solutions of the full semiclassical Einstein-Langevin equation. 
There are analogous solutions to a "reduced" version of this equation inspired in a "mini-superspace" model [ p6| , and 
there is also a previous attempt to obtain a solution to the Einstein-Langevin equation in Rcf. [ [73[ , but, there, the 
non-local terms in the Einstein-Langevin equation were neglected. 

On the other hand, the results of this calculation, which confirm our expectations that gravitational fluctuations are 
negligible at length scales larger than the Planck length, but also predict that the fluctuations are strongly suppressed 
on small scales, can be considered a first test of stochastic semiclassical gravity. In addition, we can extract conclusions 
on the possible qualitative behavior of the solutions to the Einstein-Langevin equation. Thus, it is interesting to note 
that the correlation functions at short scales are characterized by correlation lengths of the order of the Planck 
length; furthermore, such correlation lengths enter in a non-analytic way in the correlation functions. This kind of 
non-analytic behavior is actually quite common in the solutions to Langevin-type equations with dissipative terms 
and hints at the possibility that correlation functions for other solutions to the Einstein-Langevin equation are also 
non-analytic in their characteristic correlation lengths. 



A. Perturbations around Minkowski spacetime 



The Minkowski metric rj^^, in a manifold Ai which is topologically IR and the usual Minkowski vacuum, denoted 
as |0), are the class of simplest solutions to the semiclassical Einstein equation ( ^.7| ), the so called trivial solutions 
of semiclassical gravity [Q. Note that each possible value of the parameters (m^,^) leads to a different solution. In 
fact, we can always choose a renormalization scheme in which the renormalized expectation value (0| T^'^ |0)[77] — 0. 
Thus, Minkowski spacetime (H ,17^1/) and the vacuum state |0) are a solution to the semiclassical Einstein equation 
with renormalized cosmological constant A = 0. The fact that the vacuum expectation value of the renormalized 
stress-energy operator in Minkowski spacetime should vanish was originally proposed by Wald |57| and it may be 
understood as a renormalization convention [ pT|p^ . There are other possible renormalization prescriptions in which 
such vacuum expectation value is proportional to , and this would determine the value of the cosmological constant 
A in the semiclassical equation. Of course, all these renormalization schemes give physically equivalent results: the 
total effective cosmological constant, i.e., the constant of proportionality in the sum of all the terms proportional to 
the metric in the semiclassical Einstein and Einstein-Langevin equations, has to be zero. 

As we have already mentioned the vacuum |0) is an eigenstate of the total four- momentum operator in Minkowski 
spacetime, but not an eigenstate of 3^ [77]. Hence, even in the Minkowski background, there are quantum fluctuations 
in the stress-energy tensor and, as a result, the noise kernel does not vanish. This fact leads to consider the stochastic 
correction s to t his class of trivial solutions of semiclassical gravity. Since, in this case, the Wightman and Feynman 
functions (2.46), their values in the two-point coincidence limit, and the products of derivatives of two of such functions 
appearing m expressions { ^A7\ ) and ^A§i ) are known in dimensional regularization, we can compute the semiclassical 
Einstein-Langevin equation using the method outlined in section ||. 



In Minkowski spacetime, the components of the classical stress-energy tensor (2.3) reduce to 



1 



(3.1) 



where □ = 9^9'^, and the formal expression for the components of the corresponding "operator" in dimensional 
regularization, see Eq. (2.4), is 



(3.2) 



where V'^'^ is the differential operator (2.5), with g^^ — r/^^, R^^ = 0, and = d^. The field 0„(a;) is the field 
operator in the Heisenberg representation in a n-dimensional Minkowski spacetime, which satisfies the Klein-Gordon 
equation ( ^.2| ). We use here a stress-energy tensor which differs from the canonical one, which corresponds to ^ = 0, 
both tensors, however, define the same total momentum. 
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The Wightman and Feynman functions ( ^.46| ) when g^^ = rj^^, are well known: 

Gtix,y) = «A+(a; - y), Gp^{x,y) = Ap„(a; - y), 



(3.3) 



with 



A+(a;) = 



(27r)" 



(27r)« fc2 + m2-ie' 



(3.4) 



where fc^ = r/piy/c^A:" and fcx = -q^i^k^^x^ . Note that the derivatives of these functions satisfy d^A^{x—y) = 9^A+(x— y) 
and 9-^A+(a; — y) — ^d^/S.'^{x — y), and similarly fo r th e Feynman propagator Af„(x — y). 

To write down the semiclassical Einstein equation (|]^) in n-dimensions for this case, we need to compute the vacuum 
expectation value of the stress-energy operator components (3^). Since, from (^T^), we have that (0|(^^j(a;)|0) = 
iAj7„(0) = «A+(0), which is a constant (independent of x), we have simply 



T'rio)M = -*(a^9''A,j(o) = -z 



(27r)" fc2 + m2 - ie 



~\4^/ V 2/ 



(3.5) 



where the integrals in dimensional regularization have been computed i n th e standard way (see Ref. [^) and where 
r(z) is the Euler's gamma function. The semiclassical Einstein equation (2.7) in n-dimensions before renormalization 
reduces now to 



A? 



SttG, 



(3.6) 



This equation, thus, simply sets the value of the bare coupling constant Ab/Gb- Note, from (3.5), that in order to 
have (0| T^" \0)[v] = 0, t he renormalized and regularized stress-energy tensor "operator" for a scalar field in Minkowski 
spacetime, see Eq. (2^), has to be defined as 



which corresponds to a renormalization of the cosmological constant 



ri— 4 



A, 



A 



K„ + 0(n-4), 



where 



Gb G TT nin-2) 



(n-4) \47r/i2 



n-4 



V47r/i2 



+ 0(n-4). 



(3.7) 



(3.8) 



(3.9) 



being 7 the Euler's constant. In the case of a mass less scalar field, m? — Q, one simply has Ab/Gb = A/G. Introducing 
this renormalized coupling constant into Eq. (3.6), we can take the limit n^4. We find again that, for (IR^,77a6, |0)) 
to satisfy the semiclassical Einstein equation, we must take A = 0. 

We can now write down the Einstein-Langevin equations for the components /i^^ of the stochastic metric perturba- 
tion in dimensional regularization. In our case, using {0\4>l{x)\0) = iAp„{0) and the explicit expression of Eq. ( |2.43| ) 
we obtain 



1 



SnGi 



G^''"'" +Ab ( h'"'-^r|^"'h 



(x) - - aBD''''"'{x) - 2pBB^'^^''ix) 



-CG<^'''''(x)M-("-')*A,„(0)-H2yd"yM-("-4)i/r"''(2;,y)/ia/5(j;) = 2e^''(x). (3.10) 

The indices in /i^j^ are raised with the Minkowski metric and h = h^, and here a superindex (i) denotes the components 
of a tensor linearized around the flat metric. 

Note that in n-dimensions the correlator for the field is written as 
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{n^)C^{y))s = A^-^("~^)7Vr"^(x, y), (3.11) 
Explicit expressions for D'l""' and _B(i>^'' are given by 

D^'''"'{x) = I {ST^J^f - J'rK^) h^x), B^'^^^-^ix) = 2T>i'^TfKp{x). (3.12) 
where T!^" is the differential operator T!^" = 'q^'^Ox — d!^dl^- 

B. The kernels in the Minkowski background 

1. The noise and dissipation kernels 



Since the two kernels (2.45) are free of ultraviolet divergencies in the limit n — s- 4, we can deal directly with 



the F^''"'3(a; - y) = lim„^4 A^"^^""'^' in Eq. ^A§. The kernels ANt'''°'f^{x,y) = ReF^"'°''^{x - y) and 
AHj^'^"^ {x , y) — Im F'^'^"^ {x ~ y) are actually the components of the "physical" noise and dissipation kernels that will 
appear in the Einstein-Langevin equations once the renormalization procedure has been carried out. The bi-tensor 
pt^i^ap expressed in terms of the Wightman function in four spacetime dimensions, in the following way: 

F^'^^f^ix) = -2 \d''d^°'A+{x)d^^d''A+{x)+V'"'{d"'A+{x)d^A+{x)) 

+ V°'f^{d^'A+{x)d''A+{x)) +V'"'V°''^{A+^{x))] . (3.13) 
The different terms in Eq. ( 3.13| ) can be easily computed using the integrals 



I{p) ^ J ^ S{k' + m') 0{-k°) S[ik - vf + 0(fcO - /), (3.14) 

and /^'i - **'- (p) which are defined as the previous one by inserting the momenta fc^^ . . . k^'^ with r = 1, 2, 3, 4 in the 
integrand. All these integral can be expressed in terms of /(p); see Ref. [|j for the explicit expressions. It is convenient 
to separate /(p) in its even and odd parts with respect to the variables p'^ as 

/(p) = /s(p) + /a(p), (3.15) 

where /s(— p) = /s(p) and /a(— p) = — /a(p). These two functions are explicitly given by 



-1 



= signp" 0( V - 4m2) y 1 + 4 (3-16) 

After some manipulations, we find 

+ ^^r^:^/(gie-- (3AS+^)%(P), (3.17) 

where = ^—1/6. The real and imaginary parts of the last expression, which yield the noise and dissipation kernels, 
are easily recognized as the terms containing Is{p) and Ia(j>), respectively. To write them explicitly, it is useful to 
introduce the new kernels 



2887ry (27r)4'' " ^ V " ' ' p2 ' ^2 
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19207r J (27r)4 



2887r J {2tt) 



e**"" signp° 



V -4m2) Jl 



3 AC 



and we finally get 



(3.18) 



)+.Fr^,"'3^B(a;-y;m2,A0. 



(3.19) 



Notice that the noise and dissipation kernels defined in (3.18) are actually real because, for the noise kernels, only the 
cos pa; terms of the exponentials e*f^ contribute to the integrals, and, for the dissipation kernels, the only contribution 
of such exponentials comes from the isinpx terms. 



2. The kernel H'^.""'^ {x,y) 



The evaluation of the kernel components Hi^!^"^ {x,y) is a much more cumbersome task. Since these qua ntitie s 



contain divergencies in the limit ri^4, we shall compute them using dimensional regularization. Using Eq. (2. 
these components can be written in terms of the Feynman propagator (3.4) as 



where 



Kt""''^{x) = ^^l-'^''-^U2^^'^^"A^Jx)^^^^''A^Jx) + 2V^"'{^"A^„{x)^^A^Jx)) 



+ 2V"f(^d''A,^{x)d-A,^{x)) +2V^''V''^{Al^{x)) + 
+ A,„(0) {v^-'V'^^ + r^V^"") + \ V'"'V''^ (A.„(a;)n - m^A^jO)) 
Let us define the integrals 

Up) EE 



77''''a("A^„(x) a''' + 7]''^d'^^A^jx) d"^ 

6^{x) 



d"k 



1 



(27r)" (fc^ + — ie) [{k — p)^ + — ie] 
and J{^/^" '^'' {p) obtained by inserting the momenta k^^ . . . k^'' into the previous integral, together with 

d"fc 1 



-(n-4) 



(27r)" (fc2 +to2 - ie)' 



(3.20) 



(3.21) 



(3.22) 



(3.23) 



and J"on ^"^ which are also obtained by inserting momenta in the integrand. Then, the different terms in Eq. (3.21) 
can be computed. These integrals are explicitly given in Ref. 1^. It is found that /q^ = and the remaining integrals 
can be written in terms of /o„ and J„(p). It is useful to introduce the projector P^'" orthogonal to p'' as 

.2pM^^^P^p2 (3.24) 



P 



then the action of the operator JF^'' is simply written as JF^"' J d^pe^^^ f{p) = — J d^pe^^^ f{p)p^P^^'^ where f{p) is 
an arbitrary function of p^ . 

After a rather long but straightforward calculation, we get, expanding around n = A^ 
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(47r)2 



2 



3 (n-2) 
n(n-2) 



^ ('1 fti-ioi f 13)11 _ 'ptiU'pa(3\ 
-1 on ^ 3: *^a: ^ x ^ x } 



180 



^ (a^.^f-^f)" - ^r-^:'') + ^ (6oe-ii) ^r-^H '^"w 



2 1 

j^gg V'^-' a; ■'x •' X •' X ) ^ c^-j •' x •' x 



A„(a;) ^ +0(n-4), 



where k„ has been defined in (|3.9| ), and and A„(a;) are given by 



A„(x) 



(27r)" 



(3.25) 



(3.26) 
(3.27) 



where (p{p'^) = da In |1 + ^ Q!(1 — a)|. 

The imaginary part of ( 3.25 ) [which, using ( |3.2C| ), gives the kernel components Hi^J^"^ {x,y)] can be easily 

obtained multiplying this expression by —i and retaining only the real part, (p{p'^), of the function (f>{p'^). Making use 
of this result, it is easy to compute the contribution of these kernel components to the Einstein-Langevin equation. 



C. The Einstein-Langevin equations 



With the results of the previous subsections we can now write the rt- dime nsion al E instein-Langevin equation (3.10), 
previous to the renormalization. Taking into account that, from Eqs. (|3.5|) and (3.6), 



1 m'^ ^, 
SttGb 47r^ n(n—2) 



we may finally write: 
1 



G^^^^-(x) - \ aBD^^'^-^ix) - 2/3bS<^""'(.t) + 



TT) 1 

^ in-2) ^90 



1 
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28807r2 15 



— D^'^'"'{x) + - -lOAe S<^>^''(a;) 



m 
~ ~3 



(27r)" 



Jp{x-y) 



1 + 4^ Z?<^>'"'(y) + 10 3AC+^ S<^"^^(y) 



A„(a;-y) (sD^'^^'' + 5B<^>^'')(y) I + 2 / d'^y fi^'^^-^^ H^'^^ {x,y) KM + 0{n-A) 



(3.28) 



(3.29) 



Notice that the terms containing the bare cosmological constant have canceled. These equations can now be renor- 
malized, that is, we can now write the bare coupling constants as renormalized coupling constants plus some suitably 
chosen counterterms and take the limit n— >4. In order to carry out such a procedure, it is convenient to distinguish 
between massive and massless scalar fields. We shall evaluate these two cases in different subsections. 
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1. Massive field 



In the case of a scalar field with mass to 7^ 0, we can use, as we have done in Eq. (3.8) for the cosmological constant, 
the renornialized coupling constants 1/G, a and (3 as 



1 

as -- 



1 
G 



TT ^ (n-2) 
1 1 



K„ + 0(n-4), 



(47r)2 120 



K„ + 0(n-4), 



/3B=/?+:|4«„ + 0(n-4). 



(3.30) 



Note that for conformal coupling, = 0, one has 1/Gb = 1/G and /3b = /?, that is, only the coupling constant a 
and the cosmological constant need renormalization. 
Let us introduce the two new kernels 



1 



19207r2 J (27r)4 




-ITT sign p° 0(-p^-4to^) 1 / l + 4 ■ 



1 



2887r2y (27r)4 



3AC + -:t 



3 p2 



-ITT signp 



(-P^-4to^ 





(3.31) 



where ipijP') is given by the restriction to n = 4 of expression (3.26). Substituting expressions ( |330| ) into Eq. ( p9| ), 
we can now take the limit n— >4, using the fact that, for n — 4, D<^>'"'(a;) = (3/2) A<'^"'''(a:), we obtain the semiclassical 
Einstein-Langevin equations for the physical stochastic perturbations h^^ in the four-dimensional manifold M = TR^: 



1 



■ G<''^''(x) - 2 (aA<^>^'^(x) + PB^'^^"'{x)) + 



1 



-^A<^'^''(a;)+Q-10Ae^S<^'''''(a;) 



SttG" ' '"'^ ' 28807r2 

+ Jd^y[HA{x-y;mM'''^''iy) + HBix^y;m',A0B^'>^''{y)]^2e'ix), (3.32) 

where ^'^'^ are the components of a Gaussian stochastic tensor of vanishing mean value and two-point correlation 
function {^'"' {x)^"l^ {y)) s = N^"'°'f^ {x,y), given by one of the Eqs. ( |3.19D . Note that the two kernels defined in ( |3.3l| ) 
are real and can be split into an even part and an odd part with respect to the variables x^ , with the odd terms being 
the dissipation kernels Da{x\'w?') and Db{x]w? ^ A^) defined in ( 3.18| ). In spite of appearances, one can show that 
the Fourier transforms of the even parts of these kernels are finite in the limit p^ ^0 and, hence, the kernels Hx and 
H-Q are well defined distributions. 

We should mention that in Ref. [ p8[ , the same Einstein-Langevin equations were calculated using rather different 
methods. The way in which the result is written makes difficult a direct comparison with our equations ( 3.32 ). For 
instance, it is not obvious that there is some analog of the dissipation kernels. 



2. Massless field 



In this subsection, we consider the limit m— >0 of equations (3.29). The re norm alization scheme used in the previous 
subsection becomes singular in the massless limit because the expressions ( 3.30|) for as and Pb diverge when to— >0. 
Therefore, a different renormalization scheme is needed in this case. First, note that we may separate k„ in (3.9) as 



kn + \ \n.{ni? / jj?) + 0(n — 4), where 



(n-4) \47r 



1 



J'- 

2 V47r 



+ 0(n-4). 



(3.33) 
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and that, from Eq. ( 3.26| ), we have 



hm [.^(p^) + ln(m^//i^)l 



-2 + In 



(3.34) 



Hence, in the massless hmit, equations (3.29) reduce to 
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(47r)2 



D 



(27r)" 



90 



2 



(^) 



9 









(1)^1 



{y) 



-90A^25'^>'^''(y) 



+ hm 2/d"yA'-("-4)i/^A-o;3(^^^)/^^^(y)^o(^_4)^2e''.(^)_ 



These equations can be renormahzed by introducing the renormahzed couphng constants 1/G, a and /3 as 

t2 



1 _ 1 

Gb G ' 



aB — ct 



^ ^ (K„-l) + 0(n-4), Pb^P+ -^{Hn-l) + 0{n~4:). 



(47r)2 120 



327r2 



(3.35) 



(3.36) 



Thus, in the massless hmit, the Newtonian gravitational constant is not renormahz ed an d, in the conformal coupling 
case, = 0, we have again tha t Pb = P- Introducing the last expressions into Eq. (3.35), we can take the limit n— >4. 
Note that, by making m = in (3.18), the noise and dissipation kernels can be written as 



NA{x;m^ = Q) = N{x), 
£»A(a;;m2 = 0) D{x), 



NBix;m^ = 0,AO = 60A^^N{x), 
DB(a;;m2 = 0, A^ = 60 A^^^,^^.)^ 



where 



N{x) 



1 



e*P"0(-/), 



19207r J (27r)4 
It is now convenient to introduce the new kernel 



D{x) 



19207r J {2^ 



^e'f-signp" 



(3.37) 



(3.38) 



H{x-i?) 



1 



19207r2 J (27r)4 
19207r2 e-^o+ J (27r)4 



In 



lim 



In 



ITT sign p° 9{—p^) 



(3.39) 



Again, this kernel is real and can be written as the sum of an even part and an odd part in the variables x^, where 
the odd part is the dissipation kernel D{x). The Fourier transforms (3.38) and ( 3.39| ) can actually be computed and, 
thus, in this case, we have explicit expressions for the kernels in position space. For N{x) and D{x), we get (see, for 
instance, Ref. 18911 ) 



N{x) 



1 



19207r 



(a;2)2 



S\x) 



D{x) 



1 



19207r3 



sign x 



dix^) 



Six'). 



(3.40) 



where Vf denotes a distribution generated by the Hadamard finite part of a divergent integral, see Refs. [Q for the 
definition of these distributions. The expression for the kernel H{x;fi') can be found in Refs. [ pl]|9^ ] and it is given 
by 

^(^^ - 965^ d^) ^(^^0 + - '^''^ ''^^^] 

|i ^(.tO) 0(|xh A) ^ Six') + [1-7 - ln(/iA)] S\x) 



1 



lim 



9607r2 A^o+ 

See Ref. Ml for the details on how this last distribution acts on a test function. 



(3.41) 
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Finally, the semiclassical Einstein-Langevin equations for the physical stochastic perturbations /i^t^ in the massless 
case are 



SttG " ' ' ' " w ' '-^ ' 28807r2 

where the Gaussian stochastic source components have zero mean value and 



(e^''(x)r^(y)), = lim7V^''"'3(a:,y) = 



Ar(a;-y). 



(3.42) 



(3.43) 



It is interesting to consider the conformally coupled scalar field, i.e., the case A^ = 0, of particular interest because 
of its similarities with the electromagnetic field, and also because of its interest in cosmology: massive fields become 
conformally invariant when their masses are negligible compared to the spacetime curvature. We have already men- 
tioned that for a conformally coupled, field, the stochastic source tensor must be "traceless" (up to first order in 
perturbation theory around semiclassical gravity), in the sense that the stochastic variable £^ ^ = 7] ^^^^'^'^ behaves deter- 
ministically as a vanishing scalar field. This can be directly checked by noticing, from Eq. ( 3.43| ), that, when A^ = 0, 
one has {£,jiix)£,°'^{y))s = 0, since J^^ = 3n and J-'^^Tj^ = n!F°'^. The semiclassical Einstein-Langevin equations for 
this particular case (and generalized to a spatially flat Robertson- Walker background) were first obtained in Ref . , 
where the coupling constant /3 was fixed to be zero. 

Note that the expe ctati on val ue of the renormalized stress-energy tensor for a scalar field can be o btain ed by 
identification of Eqs. ( 3.32 ) and ( |3.42| ) with the components of the physical Einstein-Langevin equation ( 2.14 ). The 
explicit expressions are given in Ref. |^ The results agree with the general form found by Horowitz | |9^ , |93| | using an 
axiomatic approach and coincides with that given in Ref. pTf . The particular cases of conformal coupling, A^==0, 
and minimal coupling, A^ = — 1/6, are also in agreement with the results for this cases given in Refs. p^- p^ , p9| , p5[ 
(modulo local terms proportional to A'-'^^^^'^ and B^'-''^'^ due to different choices of the renormalization scheme). For 
the case of a massive minimally coupled scalar field, A^ = — 1/6, our result is equivalent to that of Ref. |p6|| . 



D. Correlation functions for gravitational perturbations 



In this section, we solve the semiclassical Einstein-Langevin equations ( |3.32| ) and ( |3.42| ) for the components G<^>^'' 
of the linearized Einstein tensor. In the first subsection we use these solutions to compute the corresponding two-point 
correlation functions, which give a measure of the gravitational fluctuations predicted by the stochastic semiclassical 
theory of gravity in the present case. Since the linearized Einstein tensor is invariant under gauge transformations of 
the metric perturbations, these two-point correlation functions are also gauge invariant. Once we have computed the 
two-point correlation functions for the linearized Einstein tensor, we flnd the solutions for the metric perturbations in 
the next subsection and we show how the associated two-point correlation functions can be computed. This procedure 
to solve the Einstein-Langevin equations is similar to the one used by Horowitz , see also Ref. , to analyze the 
stability of Minkowski spacetime in semiclassical gravity. 

We first note that the tensors and B'^'^'' can be written in terms of G'^'^"^ as 

^(i)M^ = - (Jc-^'^G*''^ - J';^G''^^"'), B<^'^'' = 2J^^''G<'';;:, (3.44) 

3 



where we have used that 3n = JFJ. Therefore, the Einstein-Langevin equations ( 3.32 ) and ( |3.42| ) can be seen as 
linear integro-differential stochastic equations for the components G*^"^''. These Einstein-Langevin equations can be 
written in a unified form, in both cases for m^O and for m — 0, as 

^ G''^^'''{x)-2{dA^'>^'^x)+PB'■'^^'^x)) + [d'^y[HAix-y)A^'>''''iy)+HBix-y)B''^''''{y)] ^2^'"'{x), (3.45) 



SttG 



where the new constants a and / 3, and the k ernel s Hf^ix) and Hb{x) can be identified in each case by comparison 
of this last equation with Eqs. (3.32) and ( 3.42| ). For instance, when m = 0, we have Ha{x) = H{x;^j,^) and 



Hsix) = 60 A^^ff (x; /i^). In this case, we can use the arbitrariness of the mass scale /i to eliminate one of the 
parameters a or /3. 

In order to find solutions to Eq. (|3.45| ) , it is convenient to Fourier transform them. Introducing Fourier transforms 



with the following convention f{p) = J d'^xe ^^^f{x)) for a given field f{x), one finds, from (3.44), 
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Using these relations, the Fourier transform of the Einstein-Langevin Eq. (B.45) reads 

1/3 ( 



where 



with 



F^\p{p) ^ F,{p) SU'^p.^ + F2{p)p'P^''lla.p. 



Flip) = 1 + 167rGp2 Hj^Ip) _ 2a 



HA{p) + 'iHB{p)-2a-6p 



(3.46) 

(3.47) 
(3.48) 

(3.49) 



In the Fourier transformed Einstein-Langevin Eq. (3.47), "(p), the Fourier transform of ^^"^(x), is a Gaussian 
stochastic source of zero average and 



{r{p)C^ip'))s - {2n)U\p + p')N'^'^'^^ip), 



(3.50) 



where we have introduced the Fourier transform of the noise kernel. The exphcit expression for iV'"'"^(p) is found 
from ( |3.19| ) and (|1|) to be 
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which in the massless case reduces to 



hin^7V^''"'^(p) = 9{-p^) - [p^f (3pM("p/5)-- _ pM'^p"/?) + 60Ae2(p2)2pM'>p' 



.Q/3 



(3.51) 



(3.52) 



i. Correlation functions for the linearized Einstein tensor 



In general, we can write G*^"^"^ = (G*^"^'^)s + Gf^^^ , where Gf^^^ is a solution to Eqs. (3.4£) (or, in the Fourier 
transformed version, ( 3.47 )) with zero average. The avera ges { G^^^^^)a must be a solution of the linearized semiclassical 
Einstein equations obtained by averaging Eqs. (3.45) [or (|3.47|) ]. Solutions to these equations (specially in the massless 
case, m — Q) have been studied by several authors [P2| , p3y97f |99| , |95| , p7| , particularly in connection with the problem 
of the stability of the ground state of semiclassical gravity. The two-point correlation functions for the linearized 
Einstein tensor are defined by 



e'^''"^(x,x') = (G<^"^'^(a;)G<^'"^(a;')>. - (G<^'^"(:e)), (G<^'"'3(a;')). = {G'P^"" {x)Gf^''\x')) s 



(3.53) 



Now we shall seek the family of solutions to the Einstein-Langevin equations which can be written as a linear 
functional of the stochastic source and whose Fourier transform, G*^'^'^(p), depends locally on ^"^(p). Each of such 
solutions is a Gaussian stochastic field and, thus, it can be completely characterized by the averages (G'^*'^'')^ and 
the two-point correlation functions ( |3.53 ). For such a family of solutions, Gf''^'^{p) is the most general solution to 
Eq. (3.47) which is linear, homogeneous and local in ^"''(p). It can be written as 



G<^»'''(p) = 167rGi^^''„^(p)r''(p), 



(3.54) 



where F)^'^^p(p) are the components of a Lorentz invariant tensor field distribution in Minkowski spacetime (by 
"Lorentz invariant" we mean invariant under the transformations of the orthochronous Lorentz subgroup; see Ref. [9^ ] 
for more details on the definition and properties of these tensor distributions), symmetric under the interchanges a<->/? 
and n^v, which is the most general solution of 
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(3.55) 



In addition, we must impose the conservation condition to the solutions: pj/G[^'^''(p) — 0, where this zero must be 
understood as a stochastic variable which behaves deterministically as a zero vector field. We can write D^^^p{p) = 
Dp"^ ^p{p) + D'^" ^^(p), where Dj^^afsip) is a particular solution to Eq. ( [3.55 ) and Df^ ^^{p) is the most general solution 

to the homogeneous equation. Consequently, see Eq. ( 3.54 ), we can write &l^^^{p) = G'^^"^^(p) + &^^'^''{p). To find 
the particular solution, we try an ansatz of the form 



D';\p{p) = d,{p) SU'-p^ + d^{p)p'P^-r^^p. 
Substituting this ansatz into Eqs. ( [3.55| ), it is easy to see that it solves these equations if we take 



di{p) 



Flip) 



d2{p) 



F2{P) 



Fi{p)F^{p) 



(3.56) 



(3.57) 



with 



F3(p) = fi(p) + ip'F2{p) = 1 - 48^Gp2 Hb{p) - 2P 



(3.58) 



and where the notation [ ]r means that the zeros of the denominators are regulated with appropriate prescriptions 
in such a way that di{p) and d2{p) are well defined Lorentz invariant scalar distributions. This yields a particular 
solution to the Einstein-Langevin equations: 



G^;^^''{p) = 16nGD^;\^ 



{p)e^{p) 



(3.59) 



which, since the stochastic source is conserved, satisfies the conservation condition. Note that, in the case of a 
massless scalar field, m = 0, the above solution has a functional form analogous to that of the solutions of linearized 
semiclassical gravity found in the Appendix of Ref. Notice also that, for a massless conformally coupled field, 

TO = and A^ = 0, the second term in the right hand side of Eq. (3.56) will not contribute in the correlation fimctions 



(3.53), since in this case the stochastic source is traceless. 



Next, we can work out the general form for -D^'^^^(p), which is a linear combination of terms consisting of a Lorentz 



invariant scalar distribution times one of the products (5^^, p^P'^'^rja/s, rj^^'^rja/s, fl'^'^p^Pap, P^FjJ^ and p^P'^'^ p^ Pap- 
However, taking into account that the stochastic source is conserved, we can omit some terms in D'j-^ afjip) ^^^"^ simply 
write 



G<^>^'^(p) = 16vrGI?ra,b)r''b), 



(3.60) 



where 



K\piP) = h,{p) SU^p) + h2{p)p'P^^iio.p + V^'Vap, 



(3.61) 



with hi{p), h2{p) and /i3(p) being Lorentz invariant scalar distributions. From the fact that D'h^ ^pip) niust satisfy 
the homogeneous equation corresponding to Eq. (3.55), we find that hi{p) and hz{p) have support on the set of points 
{p^} for which Fi{p) = 0, and that /i2(p) has support on the set of points {p^} for which Fi{p) = or F3(p) = 0. 



massless conformally coupled field, m - 



:0 and A^: 



Moreover, the conservation condition for G\^"^'^{p) implies t hat t he term with /i3(p) is only allowed in the case of a 



0. From (3.5C), we get 



{&^''"'{p)C^{p'))s 



{2iTf lQ^G5\p + p') Df^'^pjp) NP^^^p). 



(3.62) 



Note, from expressions ( ^.51 ) and ( |3.52| ), that the support of N'^'^°'f^{p) is on the set of point s {p^ j for w hi ch —p ^ > 



when 171 = 0, and for which —p^ — 4to^>0 when m^O. At such points, using expressions ( p9| ), (mi), ( p9| ) and 



(3.31), it is easy to see that Fi{p) and F3{p) are always different from zero, except for some particular values of 
and /3: (a) when m = 0, Af = and ;3>0; and (b) when m^O, 0<A^<(1/12) and P = {A^/32TT^)[TT/{Gm'^) + 1/36]. 

In the case (a), F3(p) = for the set of points {p''} satisfying —p^ = l/(967rG/3); in the case (b), F3{p) = 
for {p''} such that —p^ = to^/(3A^). Hence, except for the above cases (a) and (b), the intersection of the 
supports of N'^'^'^^{p) and D'h^Ap) is an empty set and, thus, the correlation function ( 3.62 ) is zero. In the 



cases (a) and (b), we can have a contribution to (3.62) coming from the term with /i2(p) in (3.61) of the form 
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^h'paiP) NP''"l^{p)=H3{p; {C})p'^P'"'N"'^Pp{p), where ^3(^1 {C}) is the most general Lorentz invariant distribution 
satisfying F^{p) H^{p; {C}) = 0, which depends on a set of arbitrary parameters represented as {C}. However, from 
( |3.5l|) , we see that N°'^''p{p) is proportional to 9{-p'^-4m'^) (1 + 4mVp^)^^^^^ (3A^ + m^/p^Y. Thus, in the case (a), 
we have N'^^Pp{p) — {) and, in the case (b), the intersection of the supports of N°'^Pp{p) and of i?3(p; {C}) is an empty 
set. Therefore, from the above analysis, we conclude that G^^^^'^ {p) gives no contribution to the correlation functions 
( p3| ), since (G<^'""'(p) ^'"^(p'll^O, and we have simply gf'""'^ {x,x') = {C-^^p" {x)G'~^^°'f^ {x')) s, where G<,^"^'^(a;) is the 
inverse Fourier transform of ( p.59| ). 

Therefore the correlation functions (3.53) can then be computed from 

(G<^''"'b) &r^{p))s = 64 (2^)6 5Hp + p') D';%,{p) D;\^{-p) NP^"^{p). (3.63) 

It is easy to see from the above analysis that the prescriptions [ ]r in the factors Dp are irrelevant in the last expr ession 
and, thus, they can be suppressed. Taking into account that Fi{—p) — F*{p), with / = 1,2,3, we get from Eqs. ( 3.56| ) 
and (|3^) 



(G^'^'-'ip) G^r^{p'))s = 64 {2nf G" ^-^^-^ 

\Fi{P)\ 



''"''{p)~^,p'P''-N"P'p{p) 



(3.64) 



This last expression is well defined as a bi-distribution and can be easily evaluated using Eq. (3.51). The final explicit 
result for the Fourier transformed correlation function for the Einstein tensor is thus 



{G'^'^'ip) &^"'^p'))s - ^ (2^)^ G^ T^^r# 

4& \Fi{p)\ 



-p'^-4m^)Jl + A — 













P J 



+ 10(^3Ae + ^^ (p^p^'^P'^f 



l-3p' 



F2{P) 



Fsip) 



(3.65) 



To obtain the correlation functions in coordinate space, Eq. ( 3.53| ), we have to take the inverse Fourier transform 
of the above result, the final result is: 



gP'^'^ix, x') = ^G^ T^'^'^P GAix -x') + ^ G^ T^'^Tf Q^{x - x'), 



(3.66) 



with 



-p^ - 4rn^) 



GAip) = 
Qnip) = ei-p'' - 4m^ 



'1 + 4^ 1 + 4^ 



3A^ 



Fi{p)\ 



\Fi{p)\ 



l-3p 



2 F2{p) 



Fsip) 



(3.67) 



and J'^""^ = Ti'^'' - Tp^'T^'^, and where Fi{p), I = 1,2,3, are given in ( ^loj) and (|^. Notice that, for a 

massless field (m = 0) , we have 

Fi(p) = l + 16^Gp2ii-(p;^2)^ 



(1 + 180AC^)i?(p; f) - 6T 



F^{p) = -^'KG 
Fsip) = 1 - 487rGp2 \60AC^H{p;fl^) ~ 2T 



(3.68) 



with /i = /i exp(19207r^Q!) and T = (3 — 60 A^^a, and where F[{p; /i^) is the Fourier transform of H{x; jjp') given in 
(p9|). 
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2. Correlation functions for the metric perturbations 



Starting from the solut ions found for the linearized Einstei n ten sor, which are characterized by the two-point 
correlation functions ( 3.66| ) [or, in terms of Fourier transforms, ( |3.65| )], we can now solve the equations for the metric 



perturbations. Working in the harmonic gauge, d^h^'^ = (this zero must be understood in a statistical sense) where 
h^ij = h^^ — {1/2)7] ^uh'^, the equations for the metric perturbations in terms of the Einstein tensor are 

□ /i^"(.T) = -2G<^'^''(a;), (3.69) 

or, in terms of Fourier transforms, p^ht^'^{p) = 2G'*^"^''(p). Similarly to the analysis of the equation for the Einstein 
tensor, we can write h^" = (h^") s+h'^'^ , where h^" is a solution to these equations with zero average, and the two-point 
correlation functions are defined by 

n^'^'^^ix, X') EE {h^^'^ix)h''^ix'))s - {h^''{x))s{h''''{x'))s - {h';\x)hf{x'))s. (3.70) 



We can now seek solutions of the Fourier transform of Eq. ( 3.69| ) of the form h'^^ {p) = 2D[p)Gf^'^'^ [p), where D{p) 



is a Lorentz invariant scalar distribution in Minkowski spacetime, which is the most general solution oip"D{p) = 1. 
Note that, since the linearized Einstein tensor is conserved, solutions of this form automatically satisfy the harmonic 
gauge condition. As in the previous subsection, we can write D{p) = [l/p^]r+Dh{p)-i where Dh(p) is the most general 

solution to the associated homogeneous equation and, correspondingly, we have h'^^ (p) = hp" (p) + {p). However, 
since Dh{p) has support on the set of points for which p^ = 0, it is easy to see from Eq. (j3]6|) (fr om the factor 



-p-^ — 4m^)) that {hf/' {p)Gf^"'^ {p')) s — and, thus, the two-point correlation functions ( ^.7(]| ) can be computed 



from {hf" {p)hf^ {p'))s = {hp'' {p)hp'^ {p'))s- From Eq. ( 3.65| ) and due to the factor 9{—p^ — Am?), it is also easy to see 



that the prescription [ ]r is irrelevant in this correlation function and we obtain 



{K''{p)hf{p'))s = j-^ (G<^>^'^(p)G<^>"''(p')>., (3.71) 



where {p) (jp')) s is given by Eq. ( 3.65|) . The right hand side of this equation is a well defined bi- 



distribution, at least for m^O (the function provides the suitable cutofi^). In the massless field case, since the noise 
kernel is obtained as the limit m ^ of the noise kernel for a massive field, it seems that the natural prescription to 
avoid the divergencies on the lightcone p^ = is a Hadamard finite part, see Refs. for its definition. Taking this 
prescription, we also get a well defined bi-distribution for the massless limit of the last expression. 
The final result for the two-point correlation function for the field h^^'^ is: 

Att 327r 
n^''"^ix,x') = —G^T^'"'^nAix-x') + —G^T^''TfnBix-x'), (3.72) 
45 9 



where Ha(p) = [l/{p^?] Qa{p) and U^ip) = [1/(^^)2] ^^(p), with gf,{p) and g^ip) given by (3.67). The two-point 



correlation functions for the metric perturbations can be easily obtained using hp^ = hf^i, — (1/2) r]fj^i,h'^. 



3. Conformally coupled field 



For a conformally coupled field, i.e., when m — Q and = 0, the previous correlation functions are greatly 
simplified and can be approximated explicitly in terms of analytic functions. The detailed results are given in Ref. 
[p|, here we outline the main features. 



When m = and = we have that ^B(a;) = and Gxip) = 0{-p^) \Fi(p)[ 
fimctions for the Einstein tensor is 



Thus the two-point correlations 



gi""^i^{x,x') = — G^jc-^"'' 



(27r)4 |l + 167rGp2H(p;^2)|2^ 



(3.73) 



P>,)/m2], seeEq. ( |3.39D . 



where H{p, ]i^) = (19207r2)-i ln[-((/ + it 

To estime this integral let us consider spacelike separated points [x — x')^ = (0,x — x'), and define y = x — x'. We 
may now formaly change the momentum variable by the dimensionless vector s^: p^ = s'^/|y|, then the previous 



integral denominator is |1 -t- 167r(Lp/|y|)2s2i/(s)|2^ where we have introduced the Planck length Lp = vG. It is clear 
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that we can consider two regimes: (a) when Lp <C |y|, and (b) when |y| ~ Lp. In case (a) the correlation funtion, 
for the 0000 component, say, will be of the order 

^0000 /,,\ 

In case (b) when the denominator has zeros a detailed calculation carried out in Ref. shows that: 

Lp 1 



^(y) 



|yp Ll\y^^ 



which indicates an exponiential decay at distances around the Planck scale. Thus short scale fluctuations are strongly 
suppressed. 

For the two-point metric correlation the results are similar. In this case we have 

^ ^ 45 "7 (2^)4 (p2)2|l + 167rGp2H(p;/22)|2 ^ ' 



The integrand has the same behavior of the correlation function of Eq. (3.73) thus matter fields tends to supress the 
short scale metric perturbations. In this case we find, as for the correlation of the Einstein tensor, that for case (a) 
above we have. 



and for case (b) we have 



-1/0000/ \ 



It is interesting to write expression (3.74) in an alternative way. If we i ntroduce the dimensionless tensor pA^^a/S _ 
3PM("p/3)i' _ piiypafj^ where P*"" is the projector defined in Eq. ( |3.24| ), to account for the effect of the operator 
we can write 

H''--f^{x,x')^—G' , ^ P ^ (3.75) 

' ' 45 i (2^)4 |l + 16^Gp2i/(p;/22)|2 ^ ' 

This expression allows a direct com paris on with the graviton propagator for linearized quantum gravity in the 
approximatio n fou nd by Tomboulis [100|. One can see that the imaginary part of the graviton propagator leads, in 
fact, to Eq. (^J5|). 



4- Discussion 



The main results of this section are the correlation functions ( p. 66 ) and (3.72). In the case of a conformal field, the 



correlation functions of the linearized Einstein tensor have been explicitly estimated. From the exponential factors 
g-ly/i-p these results for scales near the Planck length, we see that the correlation functions of the linearized 
Einstein tensor have the Planck length as the correlation length. A similar behavior is found for the correlation 
functions of the metric perturbations. Since these fluctuations are induced by the matter fluctuations we infer that 
the effect of the matter fields is to suppress the fluctuations of the metric at very small scales. On the other hand, at 
scales much larger than the Planck length the induced metric fluctuations are small compared with the free graviton 
propagator which goes like Lp/|yp, since the action for the free graviton goes like Sh ^ J d'^x Lp^hOh 

It is interesting to note, however, that these results for correlation functions are non-analytic in their characteristic 
correlation lengths. This kind of non-analytic behavior is actually quite typical of the solutions of Langevin-type 
equations with dissipative terms. An example in the context of a reduced version of the semiclassical Einstein- 
Langevin equation is given in Ref. Js^ . 

For background solutions of semiclassical gravity with other scales present apart from the Planck scales (for in- 
stance, for matter fields in a thermal state), stress-energy fluctuations may be important at larger scales. For such 
backgrounds, stochastic semiclassical gravity might predict correlation functions with characteristic correlation lengths 
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larger than the Planck scales. It seems quite plausible, nevertheless, that these correlation functions would remain 
non-analytic in their characteristic correlation lengths. This would imply that these correlation functions could not be 
obtained from a calculation involving a perturbative expansion in the characteristic correlation lengths. In particular, 
if these correlation lengths are proportional to the Planck constant h, the gravitational correlation functions could not 
be obtained from an expansion in h. Hence, stochastic semiclassical gravity might predict a behavior for gravitational 
correlation functions different from that of the analogous functions in perturbative quantum gravity |80| ] . This is not 
necessarily inconsistent with having neglected action terms of higher order in h when considering semiclassical gravity 
as an effective theory ^t\ . 



IV. FLUCTUATIONS OF ENERGY DENSITY AND VALIDITY OF SEMICLASSICAL GRAVITY 



We now turn our attention to some basic issues involving vacuum energy density fluctuations invoking only the 
simplest spacetimes, Minkowski and Casimir. Recent years saw the beginning of serious studies of the fluctuations of 
the en ergy mome ntum tensor (EMT) T^j/ of quantum fields in spacetimes with boundaries ||20|-p^ (such as Casimir 



effect [101|) |102 61], nontrivial topology (such as imaginary time thermal field theory) or nonzero curvature (such as 
the Einstein universe) j6^. A natural measure of the strength of fluctuations is x 0i the ratio of the variance Ap^ 
of fluctuations in the energy density (expectation value of the operator minus the square of the mean p taken with 
respect to some quantum state) to its mean-squared (square of the expectation value of p): 

{p') ~ {pf Ap2 



Alternatively, we can use the quantity introduced by Kuo and Ford Bit] 



(4.1) 



Assuming a positive definite variance Ap^ > 0, then < x ^ °° ^-^d < A < 1 always, with A <C 1 falling 
in the classical domain. Kuo and Ford (KF) displayed a number of quantum states (vacuum plus 2 particle state, 
squeezed vacuum and Casimir vacuum) with respect to which the expectation value of the energy momentum tensor 
(00 component) gives rise to negative local energy density. For these states A is of order unity. From this result 
they drew the implications, amongst other interesting inferences, that semiclassical gravity (SCG) ]|2^ could become 
invalid under these conditions The validity of semiclassical gravity in the face of fluctuations of quantum fields as 
source is an important issue which has caught the attention of many authors. Amongst others Phillips and Hu (PH) 
Q hold a different viewpoint on this issue from KF. This section is a summary of their investigations on this issue. 

To begin with it may not be so surprising that states which are more quantum (e.g., squeezed states) in nature 
than classical (e.g., coherent states) [103] may lead to large fluctuations in energy density comparable to the meanQ. 



Such a condition exists peacefully with the underlying spacetime at least at the low energy of today's universe. PH 
calculated the variance of fluctuations to mean-squared ratio of a quantum field for the simplest case of Minkowski 
spacetime i.e., for ordinary quantum field theory to be A = 2/5. This is a simple counter-example to the claim of 
KF, since A ~ 0(1) holds also for Minkowski space, where SCG is known to be valid at large scales. PH do not see 
sufficient ground to question the validity of SCG at energy below the Planck energy when the spacetime is depictable 
by a manifold structure, approximated locally by the Minkowski space. To them the fluctuations to mean being of 
the order unity arises from the quantum nature of the vacuum state and says little about the compatibility of the 
field source with the spacetime the quantum field lives in^ 



^This can be seen even in the ratio of expectation values of moments of the displacement operators in simple quantum harmonic 
oscillators 

■^One should draw a distinction between quantum fields in curved spacetime QFCST and semiclassical gravity: the former is 
a test field situation with quantum fields propagating in a fixed background space while in the latter both the field and the 
spacetime are determined self-consistently by solving the semiclassical Einstein equation. The cases studied in Kuo and Ford 
] pl| as well as many others are of a test-field nature, where backreaction is not considered. So KF's criterion pertains more 
to QFTCST than to SCG, where in the former the central issue is compatibility, which is a weaker condition than consistency 
in the latter. 
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PH pointed out that one should refer to a scale (of interaction or for probing accuracy) when measuring the 
validity of SCG. The conventional belief is that when reaching the Planck scale from below, QFTCST will break down 
because, amongst other things happening, graviton production at that energy will become significant so as to render 
the classical background spacetime unstable, and the mean value of quantum field taken as a source for the Einstein 
equation becomes inadequate. To address this issue as well as the issue of the spatial extent where negative energy 
density can exist, PH view it necessary to introduce a scale in the spacetime regions where quantum fields are defined 
to monitor how the mean value and the fluctuations of the energy momentum tensor change. 

In conventional field theories the stress tensor built from the product of a pair of field operators evaluated at a 
single point in the spacetime manifold is, strictly speaking, ill-defined. Point separation is a well-established method 
which suits the present concern very well, and we will discuss this method in section ^ For here we will use a 
simpler method to introduce a scale in the quantum field theory, i.e., by introducing a (spatial) smearing function 
/(x) to define smeared field operators (j>t{fx)- Using a Gaussian smearing function (with variance cr^) PH derive 
expressions for the EM bi-tensor operator, its mean and its fluctuations as functions of cr, for a massless scalar field 
in both the Minkowski and the Casimir spacetimes. The interesting result PH find is that while both the vacuum 
expectation value and the fluctuations of energy density grow as cr ^ 0, the ratio of the variance of the fluctuations 
to its mean-squared remains a constant Xd {d is the spatial dimension of spacetime) which is independent of a. The 
measure {— Xd/ixd + 1)) depends on the dimension of space and is of the order unity. It varies only slightly for 
spacetimes with boundary or nontrivial topology. For example A for Minkowski is 2/5, while for Casimir is 6/7 (cf, 
from [^). Add to this our prior result for the Einstein Universe, A = 111/112, independent of curvature, and that 
for hot flat space |104|, we see that invariably the fluctuations to mean ratio is of the order unity. 

These results allow us to address three interrelated issues in quantum field theory in curved spacetime in the 
light of fluctuations of quantum stress energy : 1) Fluctuation to mean ratio of vacuum energy density and the 
validity of semiclassical gravity. 2) The spatial extent where negative energy density can exist and its implications for 
quantum effects of worm holes, baby universes and time travel. 3) Dependence of fluctuations on intrinsic (defined 
by smearing or point-separation) and the extrinsic scale (such as the Casimir or finite temperature periodicity). 4) 
The circumstances when and how divergences appear and the meaning of regularization in point-defined field theories 
versus theories defined at separated points and/or smeard fields. This includes also the issue of the cross term. 

We begin by defining the smeared field operators and their products and construct from them the smeared energy 
density and its fluctuations. We then calculate the ratio of the fluctuations to the mean for a flat space (Minkowski 
geometry) followed by a Casimir geometry of one periodic spatial dimension. Finally we discuss the meaning of our 
finding in relation to the issues raised above. 



A. Smeared Field Operators 



Since the field operator in conventional point-defined quantum field theory is an operator-valued distribution, 
products of field operators at a point become problematic. This parallels the problem with defining the square of a 
delta function S^{x). Distributions are defined via their integral against a test function: they live in the space dual 
to the test function space. By going from the field operator (f>{x) to its integral against a test function, (f>{f) = J 4> f, 
we can now readily consider products. 

When we take the test functions to be spatial Gaussians, we are smearing the field operator over a finite spatial 
region. Physically we see smearing as representing the necessarily finite extent of an observer's probe, or the intrinsic 
limit of resolution in carrying out a measurement at a low energy (compared to Planck scale). In contrast to the 
ordinary point-defined quantum field theory, where ultraviolet divergences occur in the energy momentum tensor, 
smeared fields give no ultraviolet divergence. This is because smearing is equivalent to a regularization scheme which 
imparts an exponential suppression to the high momentum modes and restricts the contribution of the high frequency 
modes in the mode sum. 

With this in mind, we start by defining the spatially smeared field operator 

Mf^)= J .^(<,x')/x(x')dx' (4.3) 
where /x(x') is a suitably smooth function. With this, the two point operator becomes 

^ j y0(t,x')(^(t,x")/x(x')/x(x")dxdx' (4.4) 
which is now finite. In terms of the vacuum |0) (ak |0) =0, for all k) we have the usual mode expansion 
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0(^1, xi) = J dfi {ki) ^ciki Uki(ii,xi) + a^^ ul^{ti,xi) 



with 



(4.5) 



(4.6) 



where the integration measure J dfi (ki) and the normahzation constants Ny^-^ are given for a Minkowski and Casimir 
spaces by ( 4.15 ) and ( 4.21 ) respectively. 
Consider a Gaussian smearing function 



/xo(x) 



(4.7) 



with the properties / /xp (x') dx' = 1, / x' /x^ (x') dx' = xq and / |x'p /x^ (x') dx' = 2dCT^ + |xop. Using 
y Wk,(i,x)/x,(x)dx = iVfe,e-^*-^ Jlio^ J 

AT p—it w+i ki-xi— cr^ fci^ 

we get the smeared field operator 

and their derivatives 

(dt.^t,) (/xj = / [dtJih,^')) /x, (x') dx' 



,2 2 ti LJi 



(4.8) 



(4.9) 







(Vx,0t,) (/xj = J (Vx4(ii,x')) /x, (x') dx' 
= ~i J dfi{ki)ki Nk^e 



iki-xi— cr ki —itiLUi I 2iti uJi -t ^2iki-xi 



(4.10a) 



(4.10b) 



From this we can calculate the two point function of the field which make up the energy density and their correlation 
function. Letting x = (i,x) = (t2,X2) — (ti,xi), they are given by 



Ap2 it, ^;a)^^Jdfi (ki, k2) Nl iV^ (ki • ka + cc-a) 
Setting a; = we obtain the smeared vacuum energy density at one point 

p{<7)= f d^(ki)iV,>?e-^'^''=? 



2 p-2o-^ (fcJ+fe^)-ix-(ki+k2)+it(wi+W2) 



(4.11) 



(4.12) 



(4.13) 



1. Smeared- Field Energy Density and Fluctuations in Minkowski Space 
We consider a Minkowski space x R''- with d-spatial dimensions. For this space the mode density is 

d/i(k) = / k'^^^dk dVtd^i with / dVLd-i ^ —t-^- 



(4.14) 
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and the mode function normalization constant is 

Nk, = 1/^/2'^+'^ TT'^uji. (4.15) 
We introduce the angle between two momenta in phase space, 7, via 

ki • k2 = fcifc2 cos(7) — LJiiLj2 cos{j). (4.16) 
The averages of the cosine and cosine squared of this angle over a pair of unit spheres are 

d^i / dl72COs(7) = (4.17a) 

dfli [ dl]2Cos^(7)= 9 - (4.17b) 
The smeared energy density ( 4.13| ) becomes 



r(^) 



For the fluctuations of the smeared energy density operator, we evaluate ( |4.12| ) for this space and find 
A 2, ^ 1 p p /• f (l + cos(7))'A:f fcf ^, ^, 



2(2d+3)^2dJ^ 7o Js-^-Js-^-^ g2.^(fc?+fci) 

(d + i)r(^)^ 

2(3d+4)^^d^2(d+l)r(|)2 

Putting these together we obtain for the Minkowski space 

1 + d 



which has the particular values: 



12 3 1 

— for d — 1, — for d — 3, — for d — 5, — for d = 00. 

2 '5 8 3 



2. Smeared- Field in Casimir Topology 



(4.18) 



(4.19) 



Minkowski (d) — , „ , (4.20) 

1 + 3(X 



The Casimir topology is obtained from a flat space (with d spatial dimensions, i.e., xR"^ ) by imposing periodicity 
L in one of its spatial dimensions, say, z, thus endowing it with a. R^ x R'^~^ x topology. We decompose k into a 
component along the periodic dimension and call the remaining components k^ : 

k= fk_L, — j = (k_L,/n),Z = 27r/i: (4.21a) 



= ^kl+Pni'^ (4.21b) 
The normalization and momentum measure are 



//>oo P ^ 

dij{k)^ / k'^-^dk dnd-2 y (4.21c) 



Nk, = , (4.21d) 
" V2'*L7r'i-iwi 
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With this, the energy density ( 4.13| ) becomes 



/ 1 / \ 



we can write this as the sum of the two smeared Green function derivatives 

PL (a) = (Oi |((Vi(/.t) (/x))'| Oi) + (Oi (/x)) 



Ol 



where | Ol) is the Casimir vacuum. 

Since Gl(ct)^ = Gf} + G% [i = x^x^ or zz) we see how to spht the smeared energy density into a a 
term and the finite contribution: 



PL {<j)=pt + pt 



where 



div i^div I i^div 

r(^) 



= p[a) 



and 



fin y^fin , y^fin 



- (4^)(d+3)/2;d+i Z.^ ^^'^ ^^^^P 2p + d-l (2p)! r(p+|) 

With this we define the regularized energy density 

PL.reg = hm {pL (cr) ~ P {<t)) 
cr— +0 

_ d^fi3,+ir(-|) Til) 
2 (d + i)L''+ir(f + 1) 

and get the usual results: 

TT TT^ 2 TT^ 

for a = 1, — for a = 3, — for a — 5. 



6L2 ' 90L4 ' 945L6 

For the d-dimensional Casimir geometry, the fluctuations are 

,2 oo oo 



X (cos(7) ki k2 + f 111 n2 + ijJi W2) 

UJ1UJ2 

This expression can again be written in terms of products of the Green functions derivatives used above 

^""^ ^ 2{d-l) + ^^^^^'^ + GL(a),,J 

and spht into three general terms 

A 2 / \ A 2, div . A 2, cross . a 2, fin 
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The full expression for each of these terms are given in the Appendix of [g| The first term contains only the divergent 
parts of the Green functions while the last term contains only the finite parts. This is similar to the split we used for 
the smeared energy density above. What is new here is the middle term Ap^£"°^^. This comes about from the products 
of the divergent part of one Green function and the finite part of the other. That this term arises for computations of 
the energy density fluctuations is a generic feature. We will discuss in greater detail the meaning of this term later. 
The results of Appendix A in j|] give 

Api"''' = Xd {pfn' = Xd {p{a)f , (4.31a) 
Ap^'"°- = 2x.prpi", (4.31b) 



A 2,fin d {d+\) 2 . , oi N 

= 2 (^i^rog) • (4.31c) 



From this we see the divergent and cross terms can be related to the smeared energy density via 

Apl'''^ + Api'"°- = Xd{{prf + 2pi'>i"} (4.32) 

where Xd is the function that relates the fluctuations of the energy density to the mean energy density when the 
boundaries are not present, i.e., Minkowski space. This leads us to interpret these terms as due to the vacuum 
fluctuations that are always present. With this in mind, we define the regularized fluctuations of the energy density 

Api,„g = lim (ApI (/) - Xd{{pTf + 2pi'>i"}) 

= Xd.L {pL.icgf (4.33) 



where 



2 

We also define a regularized version of the dimensionless measure A 



d(d+l) , , 

Xd,L = \ ' ■ (4.34) 



^Pi.Rog + (Pi,Rog) 2 + q! + d 



and note the values: 



— for d = 1, — for d — 3, — for d = 5, 1 for d = oo. 
2 '7 16 

Following the procedures described in Appendix B of Phillips and Hu have made two plots, Fig. 1 of A(ct, L) 



and A^^jieg versus cr/L, (which we call a' here for short); and Fig. 2 of pL,Reg and ^ Ap\ versus a' . The range 

of a' is limited to < 0.4 because going any further would make the meaning of a local energy density ill-defined, as 
the smearing of the field extends to the Casimir boundary in space. (The infrared limit also carry important physical 
meaning in reference to the structure of spacetime.) 

Let us ponder on the meaning they convey. In Fig. 1, we first note that both curves are of the order imity. But 
the behavior of A (recall that the energy density fluctuations thus defined include the cross term along with the finite 
part and the state independent divergent part) is relatively insensitive to the smearing width, whereas Ai^fleg, which 
measures only the finite part of the energy density fluctuations to the mean, has more structure. In particular, it 
saturates its upper bound of 1 around a' — 0.24. Note that if one adheres to the KF criterion |6^] one would say that 
semiclassical gravity fails, but all that is happening here is that PL.Rog = while Ap\ j^^g shows no special feature. 
The real difference between these two functions is the cross term, which is responsible for their markedly different 
structure and behavior. 

In Fig. 2, the main feature to notice is that the regularized energy density crosses from negative to positive values 
at around cr' = 0.24. The negative Casimir energy density calculated in a point-wise field theory which corresponds 
to small ranges of a' is expected, and is usually taken to signify the quantum nature of the Casimir state. As 
ct' increases we are averaging the field operator over a larger region, and thus sampling the field theory from the 
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ultraviolet all the way to the infrared regi on. At large a' finite size efi'ect begins to set in. The difference and relation 
of these two effects are explained in [ 105 |: Casimir effect arises from summing up the quantum fluctuations of ALL 
modes (as altered by the boundary), with no insignificant short wavelength contributions, whereas finite size effect 
has dominant contributions from the LONGEST wavelength modes, and thus reflect the large scale behavior. As the 
smearing moves from a small scale to the far boundary of space, the behavior of the system is expected to shift from 
a Casimir-dominated to a flnite size-dominated effect. This could be the underlying reason in the crossover behavior 

of PL, Reg- 



B. Fluctuation to Mean Ratio and Spatial Extent of Negative Energy Density 

Now that we have the results we can return to the issues raised earlier. We discuss the first two here, i.e., 1) 
Fluctuations of the energy density and validity of semiclassical gravity, 2) The spatial extent where negative energy 
density can exist. We will discuss the regularization of energy density fluctuations and the issue of the cross term in 
the next subsection. 



1. Fluctuation to Mean ratio and Validity of SCG 

From these results we see that i) the fluctuations of the energy density as well as its mean both increase with 
decreasing distance (or probing scale), while ii) the ratio of the variance of the fluctuations in EMT to its mean- 
squared is of the order unity. We view the first but not the second feature as linked to the question of the validity of 
SCG. The second feature represents something quite different, pertaining more to the quantum nature of the vacuum 
state than to the validity of SCG. 

If we adopt the criterion of Kuo and Ford |^ that the variance of the fluctuation relative to the mean-squared 
(vev taken with respect to the ordinary Minkowskian vacuum) being of the order unity be an indicator of the failure 
of SCG, then all spacetimes studied above would indiscriminately fall into that category, and SCG fails wholesale, 
regardless of the scale these physical quantities are probed. This contradicts with the common expectation that SCG 
is valid at scales below Planck energy. We believe the criterion for the validity or failure of a theory should depend 
on the range or the energy probed. The findings of Phillips and Hu |^ related here seem to confirm this: Both 
the mean (the vev of EMT with respect to the Minkowski vacuum) AND the fiuctuations of EMT increase as the 
scale deceases. As one probes into an increasingly finer scale or higher energy the expectation value of EMT grows 
in value and the induced metric fluctuations become important, leading to the failure of SCG. A generic scale for 
this to 'happen is the Planck length. At such energy densities and above, particle creation from the quantum field 
vacuum would become copious and their backreaction on the background spacetime would become important [ p4| . 
Fluctuations in the quantum field EMT entails these quantum processes. The induced metric fluctuations H render 



the smooth manifold structure of spacetime inadequate, spacetime foams [ 106 including topological transitions |107] 
begin to appear and SCG no longer can provide an adequate description of these dominant processes. This picture 
first conjured by Wheeler is consistent with the common notion adopted in SCG, and we believe it is a valid one. 



2. Extent of Negative Energy Density 



It is well known that negative energy density exists in Casimir geometry, moving mirrors, black holes and worm 
holes. Proposals have also been conjured to use the negative energy density for the design of time machines [108|. 
Our results (Figures 1, 2) provide an explicit scale dependence of the regularized vacuum energy density pL,reg and 
its fluctuations A^^^eg i specifically cr/i, the ratio of the smearing length (field scale) to that of the Casimir length 
(geometry scale). For example. Fig. 2 shows that only for a/L < 0.24 is pL,reg < 0. Recall cr gives the spatial 
extent the fleld is probed or smeared. Ordinary pointwise quantum fleld theory which probes the field only at a point 
does not carry information about the spatial extent where negative energy density sustains. These r esult s have direct 
implications on wormhole physics (and time machines, if one gets really serious about these fictions [108|). If L is the 
scale characterizing the size ('throat') of the wormhole where one thinks negative energy density might prevail, and 
designers of 'time machines' wish to exploit for 'time-travel' , the results of Phillips and Hu can provide a limit on 
the size of the probe (spaceship in the case of time-travel) in ratio to L where such conditions may exist. It could also 
provide a quantum field-theoretical bound on the probability of spontaneous creation of baby universes from quantum 
field energy fluctuations. 



34 



C. Dependence of fluctuations on intrinsic and extrinsic scales 



One may ask why the fluctuations of the energy density to its mean for the many cases calculated by PH and KF 
should be the fractional numbers as they are. Is there any simple reason behind the following features observed in 
Phillips and Hu's calculations? 

a) A = 0(1) 

b) The specific numeric values of A for the different cases. 

c) For the Minkowski vacuum the ratio of the variance to the mean-squared, calculated from the coincidence limit, 
is identical to the value of the Casimir case at the same limit for spatial point separation while identical to the value 
of a hot flat space result with a temporal point-separation. 

Point a) has also been shown by earlier calculations |6l| , |6^ , and our understanding is that this is true only for 
states of quantum nature, including the vacuum and certain squeezed states, but probably not true for states of a 
more classical nature like the coherent state. We also emphasized that this result should not be used as a criterion 
for the validity of semiclassical gravity. 

For point b), we can trace back the calculation of the fluctuations (second moment) of the energy momentum 
tensor in ratio to its mean (first moment) to the integral of the term containing the inner product of two momenta 
ki • k2 summed over all participating modes. The modes contributing to this are different for different geometries, 
e.g., Minkowski versus Casimir boundary -for the Einstein universe this enters as 3j symbols - and could account for 
the difference in the numerical values of A for the different cases. 

For point c), to begin with, it is well-known that the regularization by taking the coincidence limit of a spatial versus 
a temporal point separation will give different results. The case of temporal split involves integration of three spatial 
dimensions while the case of spatial split involve integration of two remaining spatial and one temporal dimension, 
which would give different results. The calculation of fluctuations involves the second moment of the field and is in 
this regard similar to what enters into the calculation of moments of inertia for rotating objects. We suspect that the 
difference between the temporal and the spatial results is similar, to the extent this analogy holds, to the difference 
in the moment of inertia of the same object but taken with respect to different axes of rotation. 

It may be surprising that in a Minkowski calculation the result of Casimir geometry or thermal field should appear, 
as both cases involve a scale - the former in the spatial dimension and the latter in the (imaginary) temporal dimension. 
(Both cases have the same topology x S^, with the in the (imaginary) time for the former and in the space 
for the latter.) But if we note that the results for Casimir geometry or thermal field are obtained at the coincidence 
(ultraviolet) limit, when the scale (infrared) of the problem does not intercede in any major way, then the main 
components of the calculations for these two cases would be similar to the two cases of taking the coincidence limit 
in the spatial and temporal directions in Minkowski space. 

All of these cases are effectively devoid of scale as far as the pointwise field theory is concerned. As soon as we depart 
from this limit the effect of the presence of a scale shows up. The point-separated or field-smeared results for the 
Casimir calculation shows clearly that the boundary scale enters in a major way and the result for the fluctuations and 
the ratio are different from those obtained at the coincident limit. For other cases where a scale enters intrinsically in 
the problem such as that of a massive or non-conformally coupled field it would show a similar effect in these regards 
as the present cases (of Casimir and thermal field) where a periodicity condition exists (in the spatial and temporal 
directions respectively). We expect a similar strong disparity between point-coincident and point-separated cases. 
The field theory changes its nature in a fundamental and physical way when this limit is taken. This brings us to an 
even more fundamental issue made clear in this investigation, i.e., the appearance of divergences and the meaning of 
regularization in the light of a point-separated versus a point-defined quantum field theory. 



D. Regularization in the Fluctuations of EMT and the Issue of the Cross Term 

An equally weighty issue brought to light in the study of Phillips and Hu is the meaning of regularization in the 
face of EMT fluctuations. Since the point-separated or smeared field expressions of the EMT and its fluctuations 
become available one can study how they change as a function of separation or smearing scale in addition to how 
divergences arise at the coincidence limit. Whether certain cross terms containing divergences have physical meaning 
is a question raised by the recent studies of Wu and Ford . We can use these calculations to examine these issues 
and ask the broader question of what exactly regularization means and entails, where divergences arise and why they 
need to be, and not just how they ought to be treated. 

Recall the smeared energy density fluctuations for the Casimir topology has the form 

Aplia)^ Api- + Apr^^ + Apt" (4.36) 
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with 



Ap^^X^P^) =Xci{p{cT)? (4.37a) 
ApT'- = 2xdpTpt (4.37b) 
Ap|" = 2xd,L (pi")^ + terms that vanish as cr (4.37c) 

where Xd is the ratio between the fluctuations for Minkowski space and the square of the corresponding energy density: 
Ap^ = Xd(p(o'))^. Our results show that Ap|^(f7) diverges as the width a of the smearing function shrinks to zero 
with contributions from the truly divergent and the cross terms. We also note that the divergent term Ap"^'^ is state 
independent, in the sense that it is independent of L, while the cross term Ap'^™''^ is state dependent, as is the finite 
term Ap'^". 

If we want to ask about the strength of fluctuations of the energy density, the relevant quantity to study is the 
energy density correlation function H{x,y) = {p{x)p(y)) — (p(x)) {p{y)). It is finite at a; 7^ y for a linear quantum 
theory (this happens since the divergences for {p{x)p{y)) are exactly the same as the product {p{x)) (p(y))), but 
diverges as y x, corresponding to the coincident or unsmeared limit cr — > 0. 

To define a procedure for rendering Ap|^(tT) finite, one can see that there exists choices - which means ambiguities 
in the regularization scheme. Three possibilities present themselves: The first is to just drop the state independent 
Ap'^''^. This is easily seen to fail since we are left with the divergences from the cross term. The second is to neglect 
all terms that diverge as cr — > 0. This is too rash a move since Ap'^''"''" has, along with its divergent parts, ones that 
are finite in the ct — > limit. This comes about since it is of the form p^'^pf," and the negative powers of a present in 
p^'^ will cancel out against the positive powers in p|". Besides, they yield results in disagreement with earlier results 
using well-tested methods such as normal ordering in fiat space ]6l| | and zeta-function regularization in curved space 

The third choice is the one adopted by PH [g| . For the energy density, we can think of regularization as computing 
the contribution "above and beyond" the Minkowski vacuum contribution, same for regularizing the fluctuations. So 
we need to first determine for Minkowski space vacuum how the fluctuations of the energy density are related to the 
vacuum energy density Ap^ — x (p)^- This we obtained for finite smearing. For Casimir topology the sum of the 
divergent and cross terms take the form 

Api- + Apr^^ = X {(P^)' + 2pi'>f} = x{iPLf - (pi")'} (4.38) 

where x is the ratio derived for Minkowski vacuum. We take this to represent the (state dependent) vacuum con- 
tribution. What we find interesting is that to regularize the smeared energy density fiuctuations, a state dependent 
subtraction must be used. With this, just the a ^ limit of the finite part Ap|" is identified as the regularized 
fluctuations Ap^ j^gg. The ratio xl thus obtained gives exactly the same result as derived by Kuo and Ford for d = 3 
via normal ordering pl| ] and by PH for arbitrary d via the C- function . 

That thisprocedure is the one to follow can be seen by considering the problem from the point separation method, 
see section ^ For this method, the energy density expectation value is defined as the x' ^ x limit of 

p{x,x') =V,^,,Gix,x') (4.39) 

for the suitable Green function G{x, x') and 'Dx,x' is a second order differential operator. In the limit x' — )■ a;, G(a;, x') 
is divergent. The G reen function is regularized by subtracting from it a Hadamard form G^{x,x'): Giicg{x,x') — 



G(x,x') - G^{x,x') ||09[. With this, the regularized energy density can be obtained 

PRcg(a;) = lim [Vx.x'GYicgix.x')) (4.40) 

X — yx 

Or, upon re-arranging terms, define the divergent and finite pieces as 

G'^^ix, x') = G^{x, x'), G^'\x, x') = Gncg{x, x') = G{x, x') - G^{x, x') (4.41) 



and 



p{x, x') = p'^'^x, x') + p^'^ix, x') (4.42) 
p''"{x,x')=Vx^x'G'''^x,x') and pf^"(a;, x') = I?.,.'G«"(x, x') 
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so that PRcg(a;) = ^^^^x'^x p^'^{x, x'), which corresponds to the cr — s- hmit in PH's computation of the Casimir energy 
density. 

Now turning to the fluctuations, we have the point-separated expression for the correlation function 

H{x,y) =lim \\mV,,,,Vy,y,G{x,x',y,y') (4.43) 

X — >x y — ry 

where G{x,x' ,y,y') is the suitable four point function. For linear theories we can use Wick's Theorem to express 
this in terms of products of Green functions G{x, x' ,y, y') = G(x . y)G(x' ,y') + permutations of (x, x' , y, y'). Excluded 



from the permutations is G{x,x')G{y,y'). (Details are in which includes correct identifications of needed 

permutations and Green functions.) The general form is 

H{x,y) = lim \im 'Dx.x'1^y,y'G{x,y)G{x' ,y') + permutations (4.44) 

x' — >x y' — ^y ' ' 

The {x' ,y') {x,y) limits are only retained to keep track of which derivatives act on which Green functions, but 
we can see there are no divergences for y ^ x. However, to get the point-wise fluctuations of the energy density, 
the divergences from limj,^^; G{x^ y) will present a problem. Splitting the Green function into its flnite and divergent 
pieces, we can recognize terms leading to those we found for Ap|(cr): 

Hix, y) = iJ^-(x, y) + H"°^%x, y) + H^^{x, y) (4.45) 

where 

ff<^'^(x,y)= lim lim I?,,,-I?,,,.G"^'^(x,2;)G'^'^(a;',y') (4.46a) 

H''°'%x, y) = 2 lim lim V,,,,Vy^y,G'''^ {x, y)G«"(a;', y') (4.46b) 



^div ( ^ ^ ,^y^fin ( 

x' y' —'■y 



H'^'^ix.y) = \iui \imV,,,,Vy^y,G*''\x,y)G*''\x',y'), (4.46c) 

x' — >x y' — >y 

plus permutations. Thus we see the origin of both the divergent and cross terms. When the un-regularized Green 
function is used, we must get a cross term, along with the expected divergent term. If the fluctuations of the energy 
density is regularized via point separation, i.e. G{x, x') is replaced by G-R,eg{x, x') — G*'"(x, y), then we should do the 
same replacement for the fluctuations. When this is done, it is only the flnite part above that will be left and we can 
define the point-wise fluctuations as 

ApLg= limff^"(x,y) (4.47) 

S y^x 

The parallel with the smeared- field derivation presented here can be seen when the analysis of Gl{(t) and 
GL{cr) ^^^^ (given in the Appendix of Q) is considered. There it is shown they are derivatives of Green functions and 

can be separated into state-independent divergent part and state-dependent finite contribution: Gl{<j) i = G^"J-|-Gf_" , 
same as the split hereby shown for the Green function. 

When analyzing the energy density fluctuations, discarding the divergent piece is the same as subtracting from 
the Green function its divergent part. If this is done, we also no longer have the cross term, just as viewing the 
problem from the point separation method outlined above. We feel this makes it problematic to analyze the cross 
term without also including the divergent term. At the same time, regularization of the fluctuations involving the 
subtraction of state dependent terms as realized in t his c alculation raises new issues on regularization which merits 



further investigations. In a recent work Wu and Ford |11C] showed a connection between the cross term and radiation 
pressure. 



V. NOISE KERNEL, STRESS-ENERGY BI-TENSOR AND POINT SEPARATION 

As pointed out by one of us before the stress energy bi-tensor could be the starting point for a new quantum 
field theory constructed on spacetimes with extended structures. But for comparison with ordinary phenomena at low 
energy we need to find out if it behaves normally in the limit of ordinary (point-defined) quantum field theory. The 
method of point-separation is best suited for this purpose. The task is to seek a regularized noise-kernel for quantum 
fields in general curved spacetimes upon taking the coincidence limit. This was carried out by Phillips and Hu 
The following is a summary of their work. 
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PH began with a discussion of the procedures for deahng with the quantum stress tensor bi-operator at two separated 
points and the noise kernel and end with a general expression for the noise kernel in terms of the quantum field's 
Green function and its covariant derivatives up to the fourth order . (The stress tensor involves up to two covariant 
derivatives.) This result holds for x ^ y without recourse to renormalization of the Green function, showing that 
Nabc'd'{x,y) is always finite for x ^ y (and off the light cone for massless theories). In particular for a massless 
conformally coupled free scalar field on a four dimensional manifold they computed the trace of the noise kernel at 
both points and found this double trace vanishes identitically. This implies that there is no stochastic correction to 
the trace anomaly for massless conformal fields, in agreement with results arrived at in Refs. ^1 73,^. 

Now to obtain the point-defined quantities, one needs t o deal with the divergences in the coincidence limit. For 
this PH adopted the "modified" point separation scheme | 109| ,111,112 to get a regularized Green function. In this 
procedure, the naive Green function is rendered finite by assuming the divergences present for y ~* x are state 
independent and can be removed by subtraction of a Hadamard form. They showed that the noise kernel in the y ^ x 
limit is meaningful for an arbitrary curved spacetime by explicitly deriving a general expression for the noise kernel 
and its coincident form. 

After following these expositions of PH, we will end our discussion with a reflection on the issues related to regu- 
larization we started addressing in the previous section, reiterating the important role the point-separated quantities 
can play in a new approach to a quantum theory of gravity. There is a fundamental shift of viewpoint in the nature 
and meaning of the point separation scheme: in the 70's it was used as a technique (many practitioners may still view 
it as a trick, even a clumpsy one) for the purpose of identifying the ultraviolet divergences. Now in the new approach 
(as advocated by one of us PJ55[|) we want to use the point-separated expressions to construct a quantum field theory 
for extended spacetimes. 



A. Point Separation 

The point separation scheme introduced in the 60's by DeWitt |113(| was brought to more popular use in the 70's 



in the context of quantum field theory in curved spacetimes |114.5q] as a means for obtaining a finite quantum stress 
tensor^. Since the stress-energy tensor is built from the product of a pair of field operators evaluated at a single 
point, it is not well-defined. In this scheme, one introduces an artificial separation of the single point a: to a pair of 
closely separated points x and x' . The problematic terms involving field products such as ipixY becomes 4>{x)4>{x'), 
whose expectation value is well defined. If one is interested in the low energy behavior captured by the point-defined 
quantum field theory - as the effort in the 70's was directed - one takes the coincidence limit. Once the divergences 
present are identified, they may be removed (regularization) or moved (by renormalizing the coupling constants), to 
produce a well-defined, finite stress tensor at a single point. 

Thus the first order of business is the construction of the stress tensor and then derive the symmetric stress-energy 
tensor two point function, the noise kernel, in terms of the Wightman Green function. In this section we will use the 
traditional notation for index tensors in the point-separation context. 



1. n-tensors and end-point expansions 



An object like the Green function G{x, y) is an example of a bi-scalar: it transforms as scalar at both points x and 
y. We can also define a bi-tensor Ta-^...a^b'-^---b' {x,y)' upon a coordinate transformation, this transforms as a rank n 
tensor at x and a rank m tensor at y. We will extend this up to a quad-tensor Tai---a„^b'^---b'^^c'-[---c'j^^d'l' ■■■d'^'^ which 
has support at four points x, y, x' , y' , transforming as rank ni, n2, na, tensors at each of the four points. This also 
sets the notation we will use: unprimed indices referring to the tangent space constructed above a;, single primed 
indices to y, double primed to x' and triple primed to y' . For each point, there is the covariant derivative Vq at that 



^We know there are several regularization methods developed for the removal of ultraviolet divergences in the stress energy 
tensor of quantum fields in curved spacetime [po|-|2^. Their mutual relations are known, and discrepancies explained. This 
formal structure of regularization schemes for quantum fields in curved spacetime should remain intact as we apply them to the 
regularization of the noise kernel in general curved spacetimes. Specific considerations will of course enter for each method. But 
for the methods employed so far, such as zeta-function, point separation, dimensional, smeared-field |62|9| applied to simple 
cases (Casimir, Einstein, thermal fields) there is no new inconsistency or discrepancy. 
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point. Covariant derivatives at different points commute and the covariant derivative at, say, point x', does not act 
on a bi-tensor defined at, say, x and y: 

Tab';c;d' = Tab';d';c and Tab':c"^0. (5.1) 

To simplify notation, henceforth we will eliminate the semicolons after the first one for multiple covariant derivatives 
at multiple points. 

Having objects defined at different points, the coincident limit is defined as evaluation "on the diagonal", in the 
sense of the spacetime support of the function or tensor, and the usual shorthand [G{x,y)] = G{x,x) is used. This 
extends to n-tensors as 



T 



^ai---a,i. bi---bn^ ci---Cn^ di---dn.: (^■^) 



ga"' 



i.e., this becomes a rank (ni +712 + 713+ n^) tensor at x. The multi-variable chain rule relates covariant derivatives 
acting at different points, when we are interested in the coincident limit: 

[Tai---a„,b[---b'^];c = [Ta^-.-a^ b[---b'^-c] + [Tai---a,„ b[---b'^-c'] ■ (5.3) 

This result is referred to as Synge's theorem in this context; we follow FuUing's ^ discussion. 

The bi-tensor of parallel transport ga'' is defined such that when it acts on a vector Vb' at y, it parallel transports 
the vector along the geodesies connecting x and y. This allows us to add vectors and tensors defined at different 
points. We cannot directly add a vector Va at x and vector Wa' at y. But by using ga'' , we can construct the sum 

v" + ga'' Wb' ■ We will also need the obvious property 

The main bi-scalar we need is the world function a{x,y). This is defined as a half of the square of the geodesic 
distance between the points x and y. It satisfies the equation 

a = (5.4) 

Often in the literature, a covariant derivative is implied when the world function appears with indices: cr" = ct'", 
i.e. taking the covariant derivative at x, while tr" means the covariant derivative at y. This is done since the vector 
—a"" is the tangent vector to the geodesic with length equal the distance between x and y. As (T° records information 
about distance and direction for the two points this makes it ideal for constructing a series expansion of a bi-scalar. 
The end point expansion of a bi-scalar S{x, y) is of the form 

S{x, y) - + + fjf (7^42) + aVV^^), + aPfi WM^,^),, + • • • (5.5) 

where, following our convention, the expansion tensors A^ai---a,i with unprimcd indices have support at x and hence the 
name end point expansion. Only the symmetric part of these tensors contribute to the expansion. For the purposes 
of multiplying series expansions it is convenient to separate the distance dependence from the direction dependence. 
This is done by introducing the unit vector p°- = a"- j \f2a. Then the series expansion can be written 

5(x, y) - + + + ai^^^) ^ ^2^(4) ^ . . . ^5 

The expansion scalars are related to the expansion tensors via ^("^ = 2"/^Ap"?..p^p^'i • • -pP". 

The last object we need is the VanVleck-Morette determinant D{x,y), defined as D{x,y) = — det (—a. ab')- The 
related bi-scalar 

. ( ' (5.7) 

\Vgi^}J 

satisfies the equation 

(4 - (7./) - .^a'P = (5.8) 

with the boundary condition [A^/^] = 1. 



Further details on these objects and discussions of the definitions and properties are contained in ||5^] and [115|. 
There it is shown how the defining equations for a and A^/^ are used to determine the coincident limit expression for 
the various covariant derivatives of the world function ([cr;a], [o':a6], etc.) and how the defining differential equation 

for A^/^ can be used to determine the series expansion of A^/^. We show how the expansion tensors Al"La„ are 
determined in terms of the coincident limits of covariant derivatives of the bi-scalar S{x,y). (Ref. |115| details how 
point separation can be implemented on the computer to provide easy access to a wider range of applications involving 
higher derivatives of the curvature tensors. ) 
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B. Stress Energy Bi- Tensor Operator and Noise Kernel 



Even though we believe that the point-separated results are more basic in the sense that it reflects a deeper structure 
of the quantum theory of spacetime, we will nevertheless start with quantities defined at one point because they are 
more familiar in conventional quantum field theory. We will use point separation to introduce the biquantities. The 
key issue here is thus the distinction between point-defined (pt) and point-s epa rated (bi) quantities. 

For a free classical scalar field with the action Sm[g, defined in Eq. the classical stress-energy tensor is 

Tab = (1 - 2^) (j).a <l>-b + (^2^ ~ 2) ^'''^ '1^'^ 9ab + 2$ (0;p^ " Ip-abQab) 

+ 4>^ e (Rab - ^Rgabj - <\?- 9ab, (5.9) 



which is equivalent to the tensor of Eq. (2.3) but written in a slightly different form for convenience. When we make 
the transition to quantum field theory, we promote the field 0(x) to a field operator 4>{x). The fundamental problem 
of defining a quantum operator for the stress tensor is immediately visible: the field operator appears quadratically. 
Since (j){x) is an operator-valued distribution, products at a single point are not well-defined. But if the product is 
point separated (0^(x) — > 4){x)(t>{x')), they are finite and well-defined. 

Let us first seek a point-separated extension of these classical quantities and then consider the quantum field 
operators. Point separation is symmetrically extended to products of covariant derivatives of the field according to 

(-^.a) (0;6) -> \ (gj'Vp'Vb + gb^^'Va^p) (i^{x)(j){x'), (5.10) 
4> i(b;ab) \ (VaVfc + .g/ff/ Vp. V,.) ^{x)^{x'). (5.11) 

The bi-vector of parallel displacement (x,x') is included so that we may have objects that are rank 2 tensors at 

X and scalars at x' . 

To carry out point separation on (5.9), we first define the differential operator 

Tab = ^ (1 - 20 (.ga^'Va-Vb + Sb^'VaVfe.) + (^^ - ^ ffafeff"''' V, Vrf. 

-C (VaVb + ga"' gb'Va'Vb) + ^gab (v^V^ + V^'V^') 

+e (Rab - \9abR^ - ^m^gab (5.12) 

from which we obtain the classical stress tensor as 

Tab{x) = lim Tab(t>{x)4>{x'). (5.13) 

x' — >x 

That the classical tensor field no longer appears as a product of scalar fields at a single point allows a smooth 
transition to the quantum tensor field. From the viewpoint of the stress tensor, the separation of points is an artificial 
construct so when promoting the classical field to a quantum one, neither point should be favored. The product of 
field configurations is taken to be the symmetrized operator product, denoted by curly brackets: 

0(^)0(2/) - \ {k^), kv)} = \ {mky) + ky)kx)) (5.14) 

With this, the point separated stress energy tensor operator is defined as 

fab{x,x)= ^Tafc |0(x),0(x')| ■ (5.15) 

While the classical stress tensor was defined at the coincidence limit x' — > x, we cannot attach any physical meaning 
to the quantum stress tensor at one point until the issue of regularization is dealt with, which will happen in the next 
section. For now, we will maintain point separation so as to have a mathematically meaningful operator. 

The expectation value of the point-separated stress tensor can now be taken. This amounts to replacing the field 
operators by their expectation value, which is given by the Hadamard (or Schwinger) function 
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GW(x,x') = ({<^(x),0(x')}). (5.16) 
and the point-separated stress tensor is defined as 

{fab{x,x')) = ^^TabC'^Hx^x') (5.17) 

where, since Tab is a differential operator, it can be taken "outside" the expectation value. The expectation value of 
the point-separated quantum stress tensor for a free, massless (m = 0) conformally coupled (^ — 1/6) scalar field on 
a four dimension spacetime with scalar curvature R is 

{fab{x, X')) = i (gP'b G^'y,p,a + gP'a G^'l^p^b) " G^'l,p,^ gab 



12 



(5.18) 



1. Finiteness of Noise Kernel 



We now turn our attention to the noise kernel introduced in Eq. (2.11), which is the symmetrized product of the 
(mean subtracted) stress tensor operator: 



SNab,c'd-{x,y) = {[fab{x) - {fab{x))X'd'{y) - {f,^d'{y))]) 
= {\fab{x)X'd'{v)}) - 2{fab{x)){t,d'{y)) 



(5.19) 



Since Tab{x) defined at one point can be ill-behaved as it is generally divergent, one can question the soundness of 
these quantities. But as will be shown later, the noise kernel is finite for y ^ x. All field operator products present 
in the first expectation value that could be divergent are canceled by similar products in the second term. We will 
replace each of the stress tensor operators in the above expression for the noise kernel by their point separated versions, 
effectively separating the two points {x,y) into the four points {x,x\y,y'). This will allow us to express the noise 
kernel in terms of a pair of differential operators acting on a combination of four and two point functions. Wick's 
theorem will allow the four point functions to be re-expressed in terms of two point functions. From this we see that 
all possible divergences for y ^ x will cancel. When the coincidence limit is taken divergences do occur. The above 
procedure will allow us to isolate the divergences and obtain a finite result. 

Taking the point-separated quantities as more basic , one should replace each of the stress tensor operators in the 
above with the corresponding point separated version ( 5.15 ), with Tab acting at x and x' and T'd' acting at y and y' . 
In this framework the noise kernel is defined as 



8Nab,c'd'{x,y) = lim lim TabT^'d' G{x,x' ,y,y') 

x' y' — >y 



(5.20) 



where the four point function is 



G(x, x', y, y') = \ [({{^(rr), ^(x')} , [kvl ^v')}}) 

-2({0(x),0(x')})({0(y),0(y')}) 



(5.21) 



We assume the pairs {x, x') and {y, y') are each within their respective Riemann normal coordinate neighborhoods so 
as to avoid problems that possible geodesic caustics might be present. When we later turn our attention to computing 
the limit y ~* x, after issues of regularization are addressed, we will want to assume all four points are within the 
same Riemann normal coordinate neighborhood. 

Wick's theorem, for the case of free fields which we are considering, gives the simple product four point function in 
terms of a sum of products of Wightman functions (we use the shorthand notation Gxy = G-|-(x, y) = {(j){x) (f){y))): 



{(j>{x) (j){y) 4){X') 4>{y')) = Gxy' Gyx- + Gxx' Gyyl 



(5.22) 
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Expanding out the anti-commutators in (5.21) and applying Wick's theorem, the four point function becomes 
G{x,x',y,y') 

— Gxy' Gx'y Gxy Gx'y' ~\~ Gyx' Gy'x 4" Gyx Gy'x' (5.23) 

We can now easily see that the noise kernel defined via this function is indeed well defined for the limit (x', y') — > (x, y): 

G{x,x,y,y) = 2{Gly + Gi:). (5.24) 

From thi s we c an see that the noise kernel is also well defined for y ^ x] any divergence present in the first expectation 
value of (5.21) have been cancelled by those present in the pair of Green functions in the second term, in agreement 
with the results of section |l[ 



2. Explicit Form of the Noise Kernel 

We will let the points separated for a while so we can keep track of which covariant derivative acts on which 
arguments of which Wightman function. As an example (the complete calc ulati on is quite long), consider the result 
of the first set of covariant derivative operators in the differential operator ( 5.12| ), from both Tab and Tc'd' , acting on 
G{x,x' ,y,y'): 



(1 - 2^)' (g/'Vp„Wb + .g/' Vp,, Va 



9d' 



Wg>>>Wc')G{x,x',y,y') 



(5.25) 



(Our notation is that Vq acts at x, Vc' at y, Vh" at x' and Vd'" at y'). Expanding out the differential operator 
above, we can determine which derivatives act on which Wightman function: 

(1-2 0' 



'gc"'" 


9'" 


a {Gxy' ; bp'" 


Gx'y;q"d' H~ 


Gxy.bd' Gx' y' -q" 
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i.qffpi 


") 




'9' 


b {Gxy' -.ap" 


Gx' y:q" c' H 


Gxy.ac' Gx'y';q 


'p'" 








H~ Gyx' ;q 


" c' Gy' x;ap" 


' H" Gyx-ac' Gy'x 


'■.q"p' 
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(5.26) 

If we now let a;' — > a; and y' ^ y the contribution to the noise kernel is (including the factor of ^ present in the 
definition of the noise kernel): 

g ^0 {Gxy;ad' Gxy.bc' ^ Gxy-ac' Gxy-bd') 

+ (1^2^ {Gyx-ad' Gyx-bc' ~^ Gyx-ac' Gyx-bd')^ (5-27) 

That this term can be written as the sum of a part involving Gxy and one involving Gyx is a general property of the 
entire noise kernel. It thus takes the form 

Nabc'd'ix,y) = Nabc'd' [G+{x,y)]+Nabc'd' [G+{y,x)] . (5.28) 

We will present the form of the functional Nabc'd' [G] shortly. First we note, for x and y time-like separated, the above 
split of the noise kernel allows us to express it in terms of the Feynman (time ordered) Green function Gf{x, y) and 
the Dyson (anti-time ordered) Green function Goix^y): 



Nabc'd' {x,y) = Nabc'd' [GF{x,y)\ + Nabc'd' [Gnix.y)] 
The complete form of the functional Nabc'd' [G] is 



(5.29) 



*This can be connected with the zeta function approach to this problem |62| as follows: Recall when the quantum stress 
tensor fluctuations determined in the Euclidean section is analytically continued back to Lorentzian signature (r — ^ it), the 
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Nabc'd' [G] = Nabc'd' [G] + gabNc'd' [G] + gc'd'Kb [G] + gabgc'd'N [G] 



(5.31a) 



with 



8Nabc'd'[G] —(1 — 2^^ {G-c'bG-d'a + G-c'aG-d'b) + 4:£,'^ {G-c' d' G-ab + G G-abc' d' 

— 2$ (1 — 2 {G-b G-c'ad' + G-a G-c'bd' + G-^d' G-,abc' + G;c' G-abd') 
+ 2^1-20 {G-,aG-,bRc'd'+G,c'G,d' Rab) 

— 4^^ (G;afc Rc'd' + G-c'd' Rab) G + 2 Rc'd' RabG^ 



(5.31b) 



8Kb[G] =2(1-2^(^(^2^-0 G;p,6G/,+e (G;bG;p./+G;aG;p./) 

-4C(^(^2e-i^ G-P'G-abp'+i {G-^p^ G-^ab + GG-^abp'"')^ 

-(m^+^R') ((1-2 G;,G;fa-2GeG;,b) 



+2^ (^(^2^ - 2j G-yG-P + 2G^G-,p,P ) Rab 
-{m^+^R') ^RabG^ 



(5.31c) 



87V [G] = 2 ( 2e - I ) G;p.,G;P'« + 4e' (G;p.P G;g« + GG;pV 



-4^ [^^-2) [G^pG;,''"' +G-P G-,,\ 

2^-0 ((m^ + ei?) G;p- G;P' + (m^ + ^i?') G;p G'P 



-2 e ((m^ + ei?) G;p'P + (to^ + ^i?') G;^?' j G 
i (m2 + ei?) (m^ + ei?') G^ 



(5.31d) 



3. Trace of the Noise Kernel 



One of the most interesting and surprising results to come out of the investigations undertaken in the 1970's of the 



quantum stress tensor was the discovery of the trace anomaly |116|. W hen the trace of the stress tensor T = g 



is evaluated for a field configuration that satisfies the field equation (2.2) the trace is seen to vanish for massless 
conformally coupled fields. When this analysis i s car ried over to the renormalizcd expectation value of the quantum 
stress tensor, the trace no longer vanishes. Wald 1 112| showed this was due to the failure of the renormalizcd Hadamard 
function Grcn(a;,x') to be symmetric in x and x' , implying it does not necessarily satisfy the field equation (2^) in 
the variable x'. The definition of Grcn{x,x') in the context of point separation will come next.) 

With this in mind, we can now determine the noise associated with the trace. Taking the trace at both points x 
and y of the noise kernel functional ( 5.29| ): 



N[G] 



g'^'g''"' Nabc'd' [G] 



= -3GC 
9^' 



) G-.p'- 



G-7 



+^ [g.p'P' G,/ + G G-./^y ]+(^m^ + i^i?) (^ 



^ ) ^2 



time ordered product results. On the other hand, if the continuation is r —it, the anti-time ordered product results. With 
this in mind, the noise kernel is seen to be related to the quantum stress tensor fluctuations derived via the effective action as 

16iV.6c'd' = ATftc'd'L . „ ,, + Ar„\,,^,| , .,. (5.30) 
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+ ^(Rj G,p, G.P' + (^m2 + 1^^' j G,p G'f | (5.32) 
For the massless conformal case, this reduces to 

[G] = ^ {RR'G^ - 6G {RW + R'u) G + 18 ((DG) (D'G) + □'□G)} (5.33) 



which holds for any function G{x,y). For G being the Green function, it satisfies the field equation (2.2): 

aG={m^+^R)G (5.34) 

We will only assume the Green function satisfies the field equation in its first variable. Using the fact O'R — 
(because the covariant derivatives act at a different point than at which R is supported), it follows that 

□ '□G = (m^ +^i?)n'G. (5.35) 

With these results, the noise kernel trace becomes 

N[G] -^(^2 (l-30 + 3i? Q-e) ^ 

X |g2 (2™^ + ^^'^) + (i_6^) G;p'G;P' -6G^G;p'P'| 
(2 to' + R' G;p G'P - 18 e G;p G;p.PP' 

+18 Q-e) G;p>G;^''f| (5.36) 

which vanishes for the massless conformal case. We have thus shown, based solely on the definition of the point 
separated noise kernel, there is no noise associated with the trace anomaly. This result obtained in Ref. |^ is 
completely general since it is assumed that the Green function is only satisfying the field equations in its first variable; 
an alternative proof of this result was given in Ref. This condition holds not just for the classical field case, but 
also for the regularized quantum case, where one does not expect the Green function to satisfy the field equation 
in both variables. One can see this result from the simple observation used in section ||: since the trace anomaly 
is known to be locally determined and quantum state independent, whereas the noise present in the quantum field 
is non-local, it is hard to find a noise associated with it. This general result is in agreement with previous findings 
pl| , ^j73| , derived from the Feynman- Vernon influence functional formalism (T^Jt^ for some particular cases. 



C. Regularization of the Noise Kernel 



We pointed out earlier that field quantities defined at two separated points possess important information which 
could be the starting point for probes into possible extended structures of spacetime. Moving in the other (homeward) 
direction, it is of interest to see how fluctuations of energy momentum (loosely, noise) would enter in the ordinary 
(point-wise) quantum fleld theory in helping us to address a new set of issues such as a) whether the fluctuations 
to mean ratio can act as a criterion for the validity of semiclassical gravity, b) Whether the fluctuations in the 
vacuum energy density which drives inflationary cosmology violates the positive energy condition, c) How can we 
derive structure formation from quantum fluctuations, or d) General relativity as a low energy effective theory in the 
geometro- hydrodynamic limit [S^J^. 

For these inquires we need to construct regularization procedures to remove the ultraviolet divergences in the 
coincidence limit. The goal is to obtain a finite expression for t he no ise kernel in this limit. 

We can see from the point separated form of the stress tensor (5.17) what we need to regularize is the Green function 
). Once the Green function has been regularized such that it is smooth and has a well defined x' ^ x limit, 
the stress tensor will be well defined. In Minkowski space, this issue is easily resolved by a "normal ordering" 
prescription, which hinges on the existence of a unique vacuum. Unfortunately, for a general curved spacetime, there 
is no unique vacuum, and hence, no unique mode expansion on which to build a normal ordering prescription. But 
we can still ask if there is a way to determine a contribution we can subtract to yield a unique quantum stress tensor. 
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Here we follow the prescription of Wald [109|, and Adler et. al. (with corrections |112]) summarized in |2^. We 

give a short synopsis below as it will be referred to in subsequent discussions. 
The idea builds on the fact that for G{x,x')^^ — {ijj\(j3{x)(j){x')\uj) , the function 



F{x,x') = G{x,x')^,-G{x,x'\ 



(5.37) 



is a smooth function of x and x' , where uji and W2 denote two different states. This means the difference between the 
stress tensor for two states is well defined for the point separation scheme, i.e., 



1 



J^ab = 7: lim Tab {F{x, x) + F{x\ x)) 



(5.38) 



is well defined. So a bi-distribution G^{x,x') might be useful for the vacuum subtraction. At first, it would seem 
unlikely we could find such a unique bi-distribution. Wald found that with the introduction of four axioms for the 
regularized stress tensor 



(TQ6(a;))ion = lim ^-Tab (g'^^^x,x') - G^{x,x') 

x' — 'X A \ 



(5.39) 



G {x,x') is uniquely determined, up to a local conserved curvature term. The Wald axioms are |112,25 

1. The difference between the stress tensor for two states should agree with ( 5.38| ); 

2. The stress tensor should be local with respect to the state of the field; 

3. For all states, the stress tensor is conserved: V°'{Tab) — 0; 

4. In Minkowski space, the result (0|rab|0) = is recovered. 



We are still left with the problem of determining the form of such a bi-distribution. Hadamard |117| showed that 
the Green function for a large class of states takes the form (in four spacetime dimensions) 



G^{x,x')^ — 



2U{x,x') 



+ V{x, x') log a + W{x,x') 



(5.40) 



with U{x,x'), V{x,x') and W{x,x') being smooth functions^ We refer to Eqn ( 5.40D as the "Hadamard ansatz". 
Since the functions V{x,x') and W{x,x') are smooth functions, they can be expanded as 



V{x,x') 
W{x,x') = 



00 

E 

n=0 



Vn{x, x')(7' 



Wn[X, X )<7 



(5.41a) 
(5.41b) 



with the u„'s and i/;„'s themselves smooth functions. These functions and U{x,x') are determined by substituting 
G^{x, x') in the wave equation KG^{x, x') = and equating to zero the coefhcients of the explicitly appearing powers 
of (t" and cr" log a. Doing so, we get the infinite set of equations 



U{x,x')^A^/^; 
2HoVa + = 0; 

2nHnVn + KVn-1 =0, 71 > 1; 
2nHnWn + KWn~l =0, 71 > 1 



(5.42a) 
(5.42b) 
(5.42c) 
(5.42d) 



^When working in the Lorentz sector of a field theory, i.e., when the metric signature is ( — ,+,+,+), as opposed to the 
Euclidean sector with the signature (4-, -(-, -(-), we must modify the above function to account for null geodesies, since 
a{x,x') — for null separated x and x' . This problem can be overcome by using a — > cr + 2ie(t — t') + e^. Here, we will work 
only with geometries that possess Euclidean sectors and carry out our analysis with Riemannian geometries and only at the 
end continue back to the Lorentzian geometry. Nonetheless, this presents no difficulty. At any point in the analysis the above 
replacement for a can be performed. 
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with 



Hn = CT'PVp + ( n - 1 + i (Da) ) . (5.42e) 



From Eqs (5.42d), the functions Vn are completely determined. In fact, they are symmetric functions, and hence 
V{x, x') is a symmetric function of x and x' . On the other hand, the field equations only determine Wn, n> 1, leaving 
wq{x,x') completely arbitrary. This reflects the state dependence of the Hadamard form above. Moreover, even if 
wq{x, x') is chosen to be symmetric, this does not guarantee that W{x, x') will be. By using axiom (4) wo{x, x') = 0. 
With this choice, the Minkowski spacetime limit is 

where now 2a — (t — t')^ — (x — x')^ and this corresponds to the correct vacuum contribution that needs to be 
subtracted. 

With this choice though, we are left with a G'"{x,x') which is not symmetric and hence does not satisfy the field 



equation at x' , for fixed x. Wald |112] showed this in turn implies axiom (3) is not satisfied. He resolved this problem 
by adding to the regularized stress tensor a term which cancels that which breaks the conservation of the old stress 
tensor: 



{T^n = (tT) + TTTT^ffafc [vi] , (5.44) 



1 



2 



where = vi{x,x) is the coincident limit of the n = 1 solution of Eq ( 5.42d[ ). This yields a stress tensor which 



satisfies all four axioms and produces the well known trace anomaly {Ta"") — \v^\ /Stt^. We can view this redefinition 
as taking place at the level of the stress tensor operator via 

Tab Tab + 2(2r)^ g°fa ["il (^-^^^ 

Since this amounts to a constant shift of the stress tensor operator, it will have no effect on the noise kernel or 
fluctuations, as they are the variance about the mean. This is further supported by the fact that there is no noise 
associated with the trace. Since this result was derived by only assuming that the Green function satisfies th e fiel d 



equation in one of its variables, it is independent of the issue of the lack of symmetry in the Hadamard ansatz (5.40). 

Using the above formalism we now derive the coincident limit expression for the noise kernel (5.31). To get a 
meaningful result, we must work with the regularization of the Wightman function, obtained by following the same 
procedure outlined above for the Hadamard function: 

G,^n{x,y) = Grcn,+{x,y) =G+{x,y)-G^{x,y) (5.46) 

In doing this, we assume the singular structure of the Wightman function is the same as that for the Hadamard 
function. In all applications, this is indeed the case. Moreover, when we compute the coincident limit of Nabc'd', 
we will be working in the Euclidean section where there is no issue of operator ordering. For now we only consider 
spacetimes with no time dependence present in the final coincident limit result, so there is also no issue of Wick 
rotation back to a Minkowski signature. If this was the case, then care must be taken as to whether we are considering 
[Nabc'd' [G'ron,+(a;,y)]] or [Nabc'd' [G'rcn,+(y,a;)]]. 

We now have all the information we need to compute the coincident limit of the noise kernel (5.31). Since the point 
separated noise kernel Nabc'd' (x, y) involves covariant d eriv atives at the two points at which it has support, when we 
take the coincident limit we can use Synge's theorem ( ^.3| ) to move the derivatives acting at y to ones acting at x. 
Due to the long length of the noise kernel expression, we will only give an example by examining a single term. 

Consider a typical term from the noise kernel functional (5.31): 

Gycyi\c' b G-£QYi;d' a G^cyi-.c' a G^cn^d' b (5-47) 

Recall the noise kernel itself is related to the noise kernel functional via 

Nabc'd' = Nabc'd' [Grcn{x, y)] + Nabc'd' [Grcniv, x)] ■ (5.48) 

This is implemented on our typical term by adding to it the same term, but now with the roles of x and y reversed, 
so we have to consider 
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Gicnic' b Gicn:d' a ~^ ^rcn; a' d G"rcn; 6' c ~^ ^ren; c' a ^rcn ; & ~t~ G^rcn;a' c G"i-cn; 6' (5.49) 

It is to this form that we can take the coincident hmit: 

m; a' c ] [G,en;b'd] (5.50) 









a] [Gren;d'b] + 


We can now apply Syn^ 


^e's theorem: 








( [G'ren; a] 


;t/ [G^rcn;a(i]) ([Gi-cn;6_ 


;c [Gren;6c]) 




+ ([G'rcn;^ 


;a [Grcn;a(i]) ([^rcn;c 


■b — [Gren;bc]) 




H" {[Gi-cn-a_ 


;c [G'rcn;ac]) ([G^rcn;??] 


:d ~ [Grcn:bd]) 




+ ([Grcn;c_ 


;a ~ [G'rcn;ac]) ([Grcn;(i 


■ b ~ [Grcnibd]) 



(5.51) 

This is the desired form for once we have an end point expansion of Gron, it will be straightforward to compute the 



above expression. The details of such an evaluation in the context of symbolic computations can be found in [115|. 

The final result for the coincident limit of the noise kernel is broken down into a rank four and rank two tensor and 
a scalar according to 

[Nabc'd'] - Kbcd + 9abKd + 9cdN'^ + gabgcdN'". (5.52a) 
The complete expression is given in Ref. 

D. Summary Statements 

1. Further Developments 

In this section we showed how to obtain a general expression for the noise kernel, or the vacuum expectation value 
of the stress energy bi-tensor for a quantum scalar field in a general curved space time using the point separation 
method. The general form is expressed as products of covariant derivatives of the quantum field's Green function. 
It is finite when the noise kernel is evaluated for distinct pairs of points (and non-null points for a massless field). 
We also have shown the trace of the noise kernel vanishes for massless conformal fields, confirming there is no noise 
associated with the trace anomaly. This holds regardless of issues of regularization of the noise kernel. 

The noise kernel as a two point function of the stress energy tensor diverges as the pair of points are brought 
together, representing the "standard" ultraviolet divergence present in the (point-defined) quantum field theory. By 
using the modified point separation regularization method we render the field's Green function finite in the coincident 
limit. This in turn permits the derivation of the formal expression for the regularized coincident limit of the noise 
kernel. 

The general results obtained here are now applied by Phillips and Hu to compute the regularized noise kernel for 
three different groups of spacetimes: 1) ultrastatic metrics including the Einstein universe | ]6^ , hot fiat space and 
optical Schwarzschild spacetimes. 2) Robertson- Walker universe and Schwarzschild black holes, from which structure 
formation from quantum fiuctuations and backreaction of Hawking radiation on the black hole spacetimc |47| can 
be studied. 3) de Sitter and anti-de Sitter spacetimes: The former is necessary for scrutinizing primordial fiuctuations 
in the cosmic background radiation while the latter is related to black hole phase transition and AdS/CFT issues in 
string theory. 

When the Green function is available in closed analytic form, as is the case for optical metrics including hot flat 
space and Einstein Universe, conformally static spaces such as Robertson- Walker Unvierses, and maximally symmetric 



spaces such as de Sitter and Anti de Sitter , one can carry out an end point expansion according to (5.5), displaying 



the ultraviolet divergence. Subtraction of the Hadamard ansatz ( ^.40 ), expressed as a series expansion (5.41), will 



render this Green function finite in the coincident limit. With this, one can calculate the noise kernel for a variety of 
spacetimes. 

2. The Meaning of Regularization Revisited 

Thus far we have focussed on isolating and removing the divergences in the stress energy bitensor and the noise 
kernel in the coincident limit. In closing, we would like to redirect the reader's attention to the intrinsic values of 
the point-separated geometric quantities like the bitensor. As stressed earlier 0, in the stochastic gravity program. 
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point separated expression of stress energy bi-tensor have fundamental physical meaning as it contains information 
on the fluctuations and correlation of quantum fields, and by consistency with the gravity sector, can provide a 
probe into the coherent properties of quantum spacetimes. Taking this view, we may also gain a new perspective on 
ordinary quantum field theory deflned on single points: The coincidence limit depicts the low energy limit of the full 
quantum theory of matter and spacetimes. Ordinary (pointwise) quantum field theory, classical general relativity and 
semiclassical gravity are the lowest levels of approximations and should be viewed not as fundamental, but only as 
effective theories. As such, even the way how the conventional point-defined field theory emerges from the full theory 
when the two points (e.g., x and y in the noise kernel) are brought together is interesting. For example, one can ask if 
there will also be a quantum to classical transition in spacetime accompanying the coincident limit? Certain aspects 
like decoherence has been investigated before (see, e.g., |^^), but here the non-local structure of spacetime and their 
impact on quantum field theory become the central issue. (This may also be a relevant issue in noncomniutative 
geometry). The point- wise limit of field theory of course has ultraviolet divergences and requires regularization. A 
new viewpoint towards regularization evolved from this perspective of treating conventional pointwise field theory as 
an effective theory in the coincident limit of the point-separated theory of extended spacetime. 

To end this discussion, we venture one philosophical point we find resounding in these investigations reported here. 
It has to do with the meaning of a point-defined versus a point-separated field theory, the former we take as an effective 
theory coarse-grained from the latter, the point-separated theory refiecting a finer level of spacetime structure. It 
bears on the meaning of regularization, not just at the level of technical procedures, but related to finding an effective 
description and matching with physics observed at a coarser scale or lower energy. 

In particular, we feel that finding a finite energy momentum tensor (and its fiuctuations as we do here) which 
occupied the center of attention in the research of quantum field theory in curved spacetime in the 70's is only a 
small part of a much larger and richer structure of theories of fields and spacetimes . We come to understand that 
whatever regularization method one uses to get these finite parts in a point-wise field theory should not be viewed 
as universally imparting meaning beyond its specified function, i.e., to identify the divergent pieces and provide a 
prescription for their removal. We believe the extended structure of spacetime (e.g., via point-separation or smearing) 
and the field theory defined therein has its own much fuller meaning beyond just reproducing the well-recognized 
result in ordinary quantum field theory as wc take the point-wise or coincident limit. In this way of thinking, the 
divergence- causing terms are only 'bad' when they are forced to a point-wise limit, because of our present inability 
to observe or resolve otherwise . If we accord them with the full right of existence beyond this limit, and acknowledge 
that their misbehavior is really due to our own inability to cope, we will be rewarded with the discovery of new 
physical phenomena and ideas of a more intricate world. (Maybe this is just another way to appreciate the already 
well-heeded paths of string theory.) 
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FIG. 1. The dimonsionlcss fluctuation measure A = ((p'^) — (p)'^) / (p^) for the Casirnir topology, along with AL,Rcg. The 
topology is that of a flat three spatial dimension manifold with one periodic dimension of period L = 1. The smearing width a 
represents the sampling width of the energy density operator p{a). A is for the complete fluctuations, including divergent and 
cross terms, while Az,,Reg is just for the finite parts of the mean energy density and fluctuations. 



FIG. 2. The finite parts of the mean energy density p**°(c, L) and the fluctuations Ap*'°(cr, L) for the Casimir topology, as 
a function of the smearing width. 



53 




0.10 



0.20 0.30 
Smearing Width cr/L 



0.40 



Fig 1, N.G.Phillips 




Fig 2, N.G.Phillips 



